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ABSTRACT 

A finite element code is developed for analysis and 
design of three dimensional truss and frame structures. 
Structures are designed for minimum weight subject to 
constraints on: member stresses, Euler buckling, shell 

buckling, joint displacements and system natural fre- 
quencies . Structures are optimized with respect to member 
size and structure configuration. 

The finite element code may be used for analysis only, 
or may be coupled to an optimizer of the user's choice. 

The displacement method is used for static analysis, and 
structure natural frequencies are calculated via the sub- 
space iteration method. 

Flexibility is provided for expansion to other than 
tubular frame elements, and provisions are made for the 
future growth to panel and other types of structural 
elements . 

Documentation is provided to facilitate use of the 
code. A User's manual is presented with examples and re- 
sults . An explanation of how this code may be coupled to 
an optimizer is also provided. 
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INTRODUCTION 



The task of the Engineer is to provide the best solution 
to the customer's problem. Usually the "best" solution is 
the one that does the job with an adequate margin of safety, 
IS aesthetically pleasing, and is economically feasible. 

The solution may be reached through various means, but effi- 
cient use of design tools, as well as efficient use of 
materials is important since both add to the overall cost of 
the product. Design optimization is one method that can be 
used to maximize the efficiency of a structure by minimizing 
its weight and, presumably, cost. 

Optimization of structures has had continuing changes 
since its development in the early 1960 's with an active 
area of research being elastic structures. The main goal is 
to design structural systems that efficiently perform speci- 
fied purposes . Since most physical structures can be 
modeled by the finite element method, a computer program 
can be written to perform the necessary calculations to 
solve the problem. 

The ourpose of this research was to develop a finite 
element code that could be used to analyze a combined truss/ 
frame structure and could be easily coupled to an optimizer j 
thus providing a useful tool for designing such structures 



as, for example, ships' masts. This code expands the pre- 
vious work by Fitzgerald [Ref. 1] on truss structures to 
the more general six degree-of-freedom truss/frame case. 

The design problem considered in this study is the opti- 
mization of three-dimensional statically indeterminant com- 
bined truss-frame structures under varying load conditions 
and subject to a variety of constraints. The objective is 
to minimize the weight of the structure where the design 
variables are member sizes and joint coordinates . Constraints 
include: maximum normal stress; maximum joint displacements; 

minimum structure natural frequencies; Euler buckling; and 
in the case of the tubular frame elements, local or shell 
buckling . 

In the present code, all gradient information is calcu- 
lated by the finite difference method. Modification of the 
code to permit calculation of gradients analytically has 
been identified as a necessary future extension. 

This document describes the use and capabilities of the 
finite element computer code to be coupled to an optimizer. 

The user's manual presented in Chapter V contains a simple 
design example in which the program is coupled to the CONMIM 
optimization code [Ref. 2]. Additionally, guidelines for 
coupling the code to an optimizer of the user's choice are 
presented . 
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Several examples demonstrating the program under a 
variety of conditions are presented. Conclusions and recom- 
mendations for future work are given. 
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II. ANALYSIS 



A. INTRODUCTION 

When the finite element method of analysis is used to 
design optimization, two objectives must be kept in mind. 
First, the number of analyses for the structure should be 
kept to a minimum. Second, the amount of gradient informa- 
tion required during the design process should be minimized 
to shorten run times and minimize computer storage 
requirements . 

B. STATIC ANALYSIS 

Initial formulation of the problem must include approxi- 
mate member areas in the case of truss elements, and for 
frame elements, characteristic dimensions (for tubular members 
mean diameter and wall thickness); material properties (which 
may be different for each member) ; a set or sets of external 
loads; any non-structural attached masses; and specified 
joint support conditions. 

The analysis for the stresses and deflections at the 
joints must satisfy the conditions of equilibrium of forces 
at the nodes and geometric conditions of compatibility of 
deformation. In this analysis the structure is assumed to 
behave in a linearly elastic fashion. The weight of an 
individual member is not inherently included as part of the 
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specified load conditions; but, as an option, half of the 
weight of each member may be applied at the member's end- 
points as an additional load in the negative Y-direction. 

For this analysis, the following assumptions are made: 
truss and frame members are treated as discrete entities; 
truss elements have three translational degrees of freedom 
at each node and are treated as pin-connected; frame elements 
have three translational and three rotational degrees of 
freedom at each node and are treated as fixed-fixed beams; 
and loads and reactions are applied at the joints as shown 
in Figure 2.1. 

The Displacement (Stiffness) method for finite element 
analysis [Ref. 3] is utilized where 

Ku = F (3.1) 

and where K is the global stiffness matrix, F is the vector 
or vectors of applied loads , and u is the vector or vectors 
of displacements. The method used herein is an extension of 
that described by Felix and Vanderplaats [Ref. 4]. 3y 
applying the constitutive stress-displacement relationships, 
stresses in the elements may be recovered. 

C. DYNAMIC ANALYSIS 

When constraints on the system' s natural frequencies 
are to be considered, the design process requires the solu- 
tion of an eigenvalue problem. The sub-space iteration 
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Figure 2.1 FORCE ORIENTATION CONVENTION FOR .AN ELEMENT 
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method of Bathe and Wilson [Ref. 5] is used to solve for the 
desired number of lowest eigenvalues and the associated 
eigenvectors. This method is reasonably efficient for a 
small number of lowest frequencies for a large problem, and 
is well suited for re-analysis when small changes are made 
in the design. 

D. GRADIENTS 

Gradients are currently calculated with respect to member 
sizes and coordinates only by finite difference techniques . 
Inclusion of the capability to calculate gradients analytic- 
ally is identified as a needed extension to this work. 
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III. OPTIMIZATION 



A. INTRODUCTION 

The primary objective of structural optimization is to 
design systems that efficiently perform specified purposes. 
Selection of a specific optimizing algorithm must take into 
account the following: 1) the structure should be analyzed 

as few times as possible; 2) the algorithm should minimize 
the amount of gradient information required; 3) the algorithm 
should provide reasonable assurance that an optimum or near- 
optimum design will be reached. 

3. GENERAL FORMULATION 

The general statement of inequality constrained minimiza- 
tion is as follows: 



Minimize 


F(X) 




(4.1) 


Subject to: 


Gj(X) < 0 


j=l,ni 


(4.2) 


1 ^ . u 

Xi 1 X. < X. 


c 

1 — 1 
II 

•H 


(4.3) 



where F(X) is the objective function to be minimized. The 
functions Gj(X) are the set of inequality constraints to 
be met. The vector X contains the design variables. 
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The inequality constraints, Gj(X) <_ 0 .nust be 

satisfied for the design to be accepted as feasible. Side 
i u 

constraints , and , are lower and upper bounds on the 
design variables. The objective function must be minimized 
as much as possible while still satisfying all inequality 
constraints. In the case that it is not possible to satisfy 
all constraints, the most nearly feasible solution must be 
found. Felix and Vanderplaats [Ref. 4] is an excellent 
source for the basic structural design formulation. 

C. DESIGN VARIABLES 

The vector X contains the design variables; in this case, 
member characteristic dimensions for frame elements, cross- 
sectional areas for truss elements , and spacial coordinates 
of the joints. The user may elect to optimize the structure 
weight with respect to any combination of the above variables . 

For truss elements, where normal stress is dependent only 
upon the magnitude of the cross-sectional area of an element 
and not the distribution of that area, it is sufficient to 
use the area as the design variable as long as the Euler 
buckling stress can be related to the cross-sectional area. 

For frame elements where stresses are dependent on area 
distribution, more flexibility is allowed by varying two 
characteristic dimensions and allowing the code to calculate 
the section properties from these . 
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In the case of tubular elements, the characteristic 
dimensions are mean diameter and wall thickness; from which 
area, bending and polar moments, and maximum radial dimen- 
sions are generated. 

The optimum geometry problem requires that the joint 
coordinates be design variables. The x, y, and z coordinates 
of a joint are treated as separate design variables. 

In many cases it is desirable to link a set of design 
variable together to preserve symmetry, limit the number of 
variables to be solved, or to limit the number of unique 
elements to be manufactured. The code has provisions for 
design variable linking by which two or more variables may 
be linked in equality or some fixed ratio. 

D. OBJECTIVE FUNCTION 

The objective function under consideration is weight 
NE 

F(X) = I P.-A.L. (4.4) 

i=l - 1 ^ 

where is the weight density (in consistent units) of the 

material of the ith element. A,; is the cross-section area 
of the ith element, and is the length of the ith element 
and NE is the number of elements in the structure. 

E . CONSTRAINTS 

This code is designed to accommodate constraints on maxi- 
mum normal stress, Euler buckling, local or shell buckling. 
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maxinuiB joint displacements , and minim.um natural frecuencies 
or tne structure. All constraint values are normalized. 

1 . Stress : 

For truss elements , maximum normal stress is 

calculated directly from the element tensile forces . For 
frame elements; tensile stress , maximum bending and 

shear stresses and are calculated. Maximum and 

minimum normal stresses are 



0 = max magnitude of a. + ^ (^^.5) 

max ° Pi - ^ 

2 

o . = min .magnitude of a + r(a. ^ (<^.5) 

min A - A . 

n 

L 

The upper and lower bounds on stress may be different for 
each member, but are laken lo be the same for every loading 
condition . 

2 . Local or Thin Shell Buckling 

The stress ai which local or thin shell buckling 
occurs is given by: 



Gt 



0 .4E, 



n.- / t; 



(4.7) 



where the subscriot i corresoonds to the mem.ber nu.mber, F.- 



is Young's modulus 



’ m. 

"T 

t,- is the element's wall thickness. 



is the element's mean dia.meter. and 
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DESIGN VARIABLE BOUNDS 



Side constraints are imposed on the design variables 

as ; 

CdJ < CD^ _< CdV (4.3) 

U. 

CD is the characterist io dimension and CD- and CD. are 

1 1 

the minimum and maximum allowable characteristic dimensions 
of the ith element, and are taken to be the same for all 
load cases. If, as is the case in a tubular element, geo- 
metry diotates that some relationship between the design 
variables cannot be exceeded, i.e., thickness cannot exoeed 
the mean diameter, the user must arrange the bounds to pre- 
clude such an occurrence . 

G. OPTIMUM GEOMETRY 

Joint coordinates are treated as design variables with 
a seoarate design variable for each coordinate direction 
of a node. Coordinate design variables may also be linked. 
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IV . PROGRAM FEATURES 



A. INTRODUCTION 

Computer codes each have their own features and formats 
with which the user must become familiar if he is to use 
the code easily and efficiently. The SADX code developed 
in the course of this research has been designed with 
various options which are explained in this chapter. 

Chapter V contains a User's Manual with sample data for a 
typical problem that might be solved with this code: a 

truss-braced cantilever beam. This problem along with 
other numerical examples will be presented in detail with 
results in Chapter VI. 

The SADX code was written to be used as a stand alone 
analysis program or as an analysis code that could be easily 
be coupled to an optimizer (of the user's choice) through 
simple modifications to the main driver program. 

3 . FEATURES 

When the user supplies member areas , section types , 
characteristic dimensions, material data, connectivities 
and joint coordinates, along with a set of program control 
parameters; the analysis mode will calculate the weight of 
the structure. The addition of one or more sets of loading 
parameters will result in the calculation of resultant joint 
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displacements, member stresses, and/or forces for each load 
case along with the desired number of structure natural 
frequencies and m.odes . Design variables may be chosen as 
truss element areas, frame element characteristic dimen- 
sions, and joint coordinates. When coupled to an optimizer, 
the code will minimize the weight of the structure and print 
the final optimization information. The code is designed to 
be as simple to use as possible while retaining the flexi- 
bility for use on a variety of problems. 

The code's modular construction allows the user to use 
frame element cross-sections other than tubular elements. 
This is done by reading in two characteristic dimensions for 
each frame element along with a section type identifier. 
Subroutine SADX35 is called to calculate the section proper- 
ties, area, maximum radial dimensions, and bending and polar 
moments of inertia. The user may augment this subroutine to 
calculate section properties for whatever section type he 
may choose to work with. 

Many computer codes require the definition of an auxil- 
iary node to orient the principle axis of a non-axially 
symmetric element . In this code it is assumed that the 
element's local z-axis is parallel to the global x-z plane. 

Pseudo-dynamic storage is used to allow storage of most 
data in one integer and one real array for more efficient 
use of storage. 
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Various print control options are available to tailor 
the printed output to match the user's desires. 

Design variable linking is available to allow elements 
and joints to be grouped to maintain symmetry, limit the 
number of independent design variables, or reduce the 
variety of member sizes generated. 

Optimization may be performed with respect to member 
size, with respect to structure geometry, or with respect to 
both . 

Structures may be comprised of truss elements, frame 
elements, or a mix of the two types. 

Structures may be optimized for multiple load cases with 
constraints imposed upon any combination of maximum normal 
stress, maximum joint displacements and rotations, Euler 
and local buckling, and a specified number of minimum natural 
frequencies in free vibration. Separate displacement con- 
straints may be imposed for each load case . Both consis- 
tent and lumped mass options are coded. Either forces or 
stresses or both can be output. The user may decide whether 
to include the structure's weight and the fixed masses as 
loads applied to the structure. 

C . EXAMPLE 

The following example of the truss-braced cantilever 
beam presented in Tables (I-V) demonstrates some of the 
options available in the code. 



25 



TABLE I 



TRUSS -BRACED CANTILEVER BEAM: INPUT CONTROL PARAMETERS 



INPUT PARANETERS FOR STRUCTURAL ANALYSIS AND DESIGN 
ROUTINE, "SADX" 

TRUSS-ERACED CANTILEVER BEAM: EXAMPLE 1 



&NTS , 
EMENTS, 



CONTROL PARAMETERS 

TOTAL NUMBER OF ELEMENTS 
TOTAL NUMBER OF BAR ELSM 
TOTAL NUMBER OF FRAME EL 
TOTAL NUMBER OF JOINTS, 

JOINT MAXIMUM DEGREES OF 
JOINT CONSTRAINT VARIABL 
NUMBER OF MATERIAL TYPES 
NUMBER OF LOAD CONDITION 
NUMBER OF EIGENVALUES, 

NO. OF EIGENVALUES CALCULATED 
NUMBER OF FIXED MASSES, 

EULER BUCKLING CONSTRAIN 



FREEDOM, 



9 

S, 



LOCAL BUCKLING CONSTRAIN 
NO. OF DISPL. CONSTRAINT 
NO. OF FREQ. CONSTRAINTS 
LUMPED MASS OPTION, 
FORCE/STRESS PRINT OPTIO 
OPTIMUM SIZS/BOTH/GEOM 0 
STRUCTURE WEIGHT AS LOAD 
FIXED MASSES AS LOADS OP 



T 

T 

S, 



ID, 

ID, 



pillON 
S OPTION 
TION 



HE: 
NEB= 
NEF = 
NJ= 
JCM = 
NCJ= 
NMT = 
NLC = 
NEIG = 
NEIG1 = 
NFM ASS: 
NEU3C: 
LBUCK = 
ND3PLC: 
NFREQ: 
LMASS = 
NSIRES: 
IDVCLC: 
NSTW: 
NFMW: 



4 
2 
2 

5 

6 
3 
2 
2 
1 
2 
1 
1 
1 
2 
1 
0 
2 
2 
1 
1 



ACCELERATION DUE TO GRAVITY, GRAV = 0.38640 

EIGENVALUE CONVERGENCE TOLERANCE, SP3EIG = 0.10000 
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TABLE II 



EXAMPLE 1: JOINT COORDINATE AND MEMBER INPUT 



JCINT CCCROIN^TES 



JCINT 

1 

2 

3 

4 

5 



0.0 

0.1C0Q0E^Q3 

0.2CCC0E+C3 

0*0 

Q.O 



0.0 
0.0 
0* 0 

C* 15000E+03 

0.0 



CCOROIN^TE CESIGN VARIABLES 
CESIGN 
VARIABLE 
X Y 2 
C 
C 
C 

c 



JCINT 

1 

2 

3 

4 

5 



C 

C 

c 

c 

c 



0 

0 

0 

1 

0 



X 
0*0 
0*0 
0.0 

O.lCCOE+01 

C.lOCOE+01 



0.0 

0.0 

0.0 

0.0 

0.500 00S+02 



MULTIPLIER 

Y 

0.0 
0.0 
0.0 

O.lOOOE+01 
Q.lOOCE^-01 



z 

0.0 

0.0 

0.0 

0. 10006+01 
0.1 0006+01 



ELEMENT INFOPMATICN FCP EAR ELEMENTS 

ELEMENT-JOINT PEL ATIONSHI PS 

LNG NODEl NOOE2 MATL CVARl AREA 
12411 C. 3000E+C1 
22512 0.3000E+C1 



LENGTH 
0. 1303E+03 
0. 1118E+03 



ELEMENT INFORMATICN FCP FRAME ELEMENTS 
ELEMENT-JOINT RELATICNSHI PS ELEMENT PROPERTIES 
ELEMENT NUMBER 3 

NCDEl NG0E2 MATL 1ST OESVAR 2ND DESVAR 
1 2 2 3 4 

AREA LENGTH I-HOMENT Y-MOMENT 

0.157E+C2 O.lOOE+03 0.511E+02 0. 5116+02 



ZMAX 

0.300E+01 



ELEMENT NUM06F 4 

NCDEl N0DE2 MATL 1ST OESVAR 2ND CESVAR 
2 3 2 3 5 

AREA LENGTH 2-MOMENT Y-M0M6NT 

0.157E+C2 0.1CCE+C3 C.511E+02 0.511E+02 



ZMAX 

0. 3006+01 



DATA 



YMAX 

•300EH 



YMAX 
• 300EU 
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TABLE III 



EXAMPLE 1; SAMPLE DISPLACEMENT AND FORCE/STRESS OUTPUT 



. c# 


EAR ELIf>EN7$ 


• 2 


INC 


LOAO CCNO 


7ENSUE FORCE IN X 


1 


1 


-C.73C7F*OA 


1 


2 




LNC 

! 


LOAO CCKO 
1 
2 


7ENSIL£ FORCE IN * 
0.7A28I ♦OR 
-0.2A13E*05 


. OF 


FMhC EiEMEftlS • 2 






ElENEAT 


LOAD CONO 


NCOS 


TENSILE FCPCE 
IN t 01 REX 


SHEAR FORCE 
IN V OIREX 


SHEAR FORCE 
IN Z OIREX 


7 CPSICN ROMENT 
xacuT X AMIS 


8 EN 01 NC MCMENT RE»rt)IFC MCHf 
A 8 CUT Y AXIS A 8 CUT Z AXIS 


NT **i« 

1 STRESS 


3 

3 


{ 


LCh 

HIGH 


- 0 . 17 S 7 E» 0 A 

O.irsTE^OA 


O.lOAEE^OA 

- 0 . 1066 E^ 0 A 


- 0 , 339 SE* 0 l 

0 . 3395 E* 0 l 


0.0 

0.0 


0 . 3395 E^C 3 

- 0 . 97 fc 6 E -03 


§:?UU:3t 




3 

3 


2 

2 


LCV 

HIGH 


0 . 66 A 7 =^ 0 A 

- 0 . 66 A 7 E» 0 A 


- 0 .A 16 7 E 403 
0 ,AU 7£»03 


- 0 . 1 A 82 E+ 0 A 

0 «lA 82 c« 0 A 


0.0 

0.0 


0 .A 816 EAC 5 

O.LOCOE^CA 




-3:n‘4l:35 


EtEHENT 


LOAC CCNO 


KCCE 


TENSILE FCRCf 
IN X OIREX 


shear FORCE 
IN T OIREX 


SHEAR FORCE 
IN 1 OIREX 


TOKSICN FOMENT 
A 8 CUT X AXIS 


SEKOINO MOMENT 8 EN 01 XC MCMENT '»£« 

A 8 CUT Y AXIS ABOUT Z AXIS STRESS 


A 


1 

1 


LCA 

HIGH 


- 0 ,A« 8 ?EH 33 

0 .AS 82 E -03 


-0 .TIASE^OS 
0 . 71 ASE «03 


- 0 .A 768 E -05 

0 .A 768 E- 0 S 


0.0 

0.0 


0 . 1373 E -03 

0 . 732 AE -03 


- 0 . 7145 E^ 0 S 

0 .A 2 S 0 E -01 


H 3 . 2 « 72 E *08 

0 . 2271 E -01 


A 2 LCh 0 . 3906 E -02 

A 2 HIGH • 0 . 390 A £-02 

EIGENVALUES ARC EIGENVECTORS 


0 . 26 S 5 E ♦OS 
- 0 . 289 SE ^03 


C.IOOOE^OA 

-o.iooal^OA 


u 


-C.lOCO|*OA 


g. 2 eS 5 E *08 
- 0 . A 25 CE -01 


0 . 3600 F« 0 S 

- 0 . 9915^-01 



HCOE Mi'aER 
ffUUI^CY • 

ETCENVAIUE • 
ElCENVECTOft 
JCINT 1 



I0«EE»O2 C?5 
fc3lSE»03 UPS 
927CE*C5 
ftEE CF FREECOX 
2 



c;ir7ni:8i 

0.17155E-01 



0.100 00E*01 



lOAC CC^O^TIO^ 

otsplacpments 

JNT >-CISPL 

1 C.C 

2 -C.1C56-03 

3 -O.I05E-O3 

i o.c 

3 C,C 

ICAC CC^0ITIC^ 

CISPLACfHENTS 
JAT X-CISPU 
1 C.C 
^ *0*i26£-03 

3 -C.226E-02 

i 

5 0.0 



C|GPEE OP FREEOCH 



Y-CISPL 
O.C 

•0.275E-02 
•O^ 1C5E-01 
O.C 
C.C 



2-01 SPL 

0.0 

-0.1376-03 
•0. 3436-03 
0.0 
0.0 



06GR6E CP FREEDOM 



Y-CISFL 
O.C 

-0.2716-02 
•0. 1616-01 
C.C 
0.0 



2-OISPL 

0.0 

-0. 187E-02 
-0.7846-02 
0.0 
0.0 



FCT A8T X 
0.0 
C.O 
C.O 
0.0 
0.0 



ROT A8T X 

0.0 



0.0 

0.0 

0.0 

Q.O 



RCT ABT Y 

0.0 

0.2066-05 
0. 2066-05 
9.0 
0.0 



ROT ABT Y 

0.0 

0.^166-04 

0.6886-04 

0.0 

0.0 



ROT A8T Z 

0.0 

-0.5676-04 

-0.8756-04 

0.0 

0.0 



ROT ABT 2 

0.0 

- 0.847 6-04 
-0.1426-03 
0.0 
0.0 



T ABT Z 

• 0 

• 566466—04 
•874816-04 

• 0 
• 0 



T ABT 2 

• 0 

.846096-04 
. 14 2696-03 

• 0 
• 0 
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TABLE IV 



EXAMPLE 1: SAMPLE CONSTRAINT OUTPUT INFORMATION 



DISPLACEMENT CCNSTRAINTS START AT NUMBER 2 AND STOP AT NUMBER 5 

TRUSS ELEMENT STRESS CCNSTRAINTS 

START AT NUMBER 6 AND STOP AT NUMBER 17 

FRAME ELEMENT STRESS CCNSTRAINTS 

START AT NUMBER 18 AND STOP AT NUMBER 41 

TFE NUMBER CF CCNSTRAINTS NCTOT = 41 



1) -.4CE882E+C2 
5> -.174111E4C1 
R) -.101236E4C1 
13) -.S75722E+CC 
17) -.901716E4CC 
21) -.471251E+C0 
5) -.1AC2C1E+C1 
9) -.<35S641E4C0 
3) -,«;cS9«39E4C0 
7) -.1S9994S+C1 
41) -.951795E+C0 



-.161055E+O1 
-.«fil588E+00 
SS7t4CE+00 
102139E+01 
126755E+C1 
-.118559E+01 
-.22SC14E+CO 
171444E401 
123411 £+01 
S99S98E+00 



-.389445E+00 

101841E+C1 

-.103769E+01 

-.888444E+00 

-.732450E+CC 

-.972695E+00 

177499E+C1 

285557E+C0 

-.965558E+00 

iOCOOCE+Cl 



-.258895E+0C 

-.943857E+00 

-.102428E+C1 

-.111156E+01 

-.152875E+C1 

-.597988C+00 

-.127432E+01 

lOOOOOE+01 

-.581741E-04 

-.132766E+01 
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TABLE V 



EXAMPLE 1: FINAL OPTIMIZATION INFORMATION 



1 



FIN4L OPTIM2ATICN INFCRMATION 

Cej = C*5675C76+C3 

DECISION V/RIAELES (X-VECTGR) 

1) C*3^ECJE+01 O.a^S^ZE^-Ol 0*4a389E^01 
7) C.45^3iE>02 

CCNSTPAIM \ALLES (G-VECTORI 
l)-C*^03866+02 -0.161C6E+01 -0. 389A5E>0C 
7 )-C.10ia«e+01 -0.94336E+00 -0.10124E+01 
13)-C*97572c*»-00 -0.10214E+01 -0*g8844£+CG 
19)-C.73245E^-00 -0. 15237OC1 -0.47125E+CC 
25 »-C.1402CE>01 -0.22501E+00 -0.17750E401 
31 )-0. 2855^ = ^00 -C.lOOCOE^Cl -C. 100006+01 
37)-C.199<5EE+01 -O.lOOOOE+01 -O.lOOOOE+01 

THERE ARE 1 ACTIVE CONSTRAINTS 
CONSTRAINT NUMBERS ARE 
36 



0.21898E+00 



-0.25889E+00 

-0.98764E+00 

-0.111165+01 

-0.11856E+C1 

-0.12743E+01 

-0.123415+01 

-0.13277E+C1 



0.15576E+00 



-0.17411E+01 

-0.10377E+C1 

-0.901725^00 

-0.9727QE+CC 

-0.95964E+00 

-O.Q6556E+00 

-0.95160E+C0 



THERE ARE 0 VICLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

TEPMNATICN CRITERION 

AeS(CBJ( n-CBJ( I-IJ) LESS THAN OABFUN FCR 3 ITERATIONS 
NUMBER OF ITE«ATICNS * 12 

OBJECTIVE FUNCTION WAS EVALUATED 114 TIMES 

CCNSTRAINT FUNCTICNS WERE EVALUATED 114 TIMES 

THIS PUN RECUIRED 116 STRUCTURAL ANALYSES 

NUMBER OF SECCNCS REQUIRED FCR EXECUTION = 2.79 

WEIGHT OF structure GIVEN AREAS S LENGTHS 
WEIGHT^ 0.5E751E+C3 

TCTAL WEIGHT INCLUDING FIXED MASSES 
TCTAL WEIGHT* C. 63751E+03 



JCINT CQCFDINATES 
JCINT X 

1 0 .0 

2 G.1C0C0E+C3 

3 0.2CCC0E+C3 

4 0.0 

5 O.C 



0.0 

0.0 

0.0 

0.84259E+02 

0.0 



0.0 
0 .0 
0.0 
0 .0 

0.454326+02 



ELEMENT INFCRMATICN FCR 3AR ELEMENTS 
ELEMENT- JCINT R ELATI ONSHI PS 



ELEMENT 

1 

2 



NODE 1 
2 
2 



NOCE 2 

4 

5 



AREA 

C.3481E+C1 

0.8499E+01 



LENGTH 

1308E+03 

1098E+03 



ELEMENT INFCRMATICN FOR FRAME ELEMENTS 

ELEMENT-JCINT RELATIONSHIPS 

LNG NCCEl NCCE2 AREA LENGTH CHAR.OIM.l 

212 0.3329E+01 O.lOOOE+03 0.4839E+01 

423 0.2266E+01 0.1000E+C3 0.48396+01 



CHAR. DIM. 2 
0.2190E+00 
0.15586+00 



30 



V. USER GUIDE 



A. INTRODUCTION 

In developing any computer code for engineering analysis, 
it is necessary to additionally develop concise, easily 
understood documentation. This SADX USER'S GUIDE is written 
to be easily understood by the user having only minimal know- 
ledge of the FORTRAN language. The format follows that of 
the optimization code, COPES /CONHIN [Ref. 2]. 

This chapter is devoted to acquainting the user with the 
code and necessary input data. 

3. DESIGN EXAMPLE 

The simple example of a four-element combined truss /frame 
structure is used to demonstrate some of the features of the 
SADX program. 

The structure is shown in Figure 5.1, and consis's of 
two tubular frame elements along the x-axis with two truss 
braces to the y and z axes from rhe joint between the frame 
elements. A non-structural fixed mass is attached at the 
outboard end of the second frame element where two loads are 
applied . 

C . SADX DATA 

The SADX program reads data from uni” 5 and wrires out- 
put on unit 6. Units 30 and 40 are used as scratch files. 
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ORIGINAL 
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The scratch file numbers may be changed by changing two cards 
at the beginning of subroutine SADXOl. The SADX program has 
the capability to read unformatted data. The following sec- 
tions define the data which is required by SADX. The data 
is segmented into "BLOCKS" for convenience. 

SADX data begins with a TITLE card and ends with a END 
card. Comment cards may be inserted anywhere in the SADX 
data stack prior to the END card, and are identified by a 
dollar sign C$) in column 1. Data may be of either the 
"110" or "FIO.O" type or may be free-format separated by 
commas with no imbedded blanks . Lines of formatted and 
unformatted data may be intermixed. 

1 . Formatted Data 

Formats are of "110" and "FIO.O" type. "I" formats 
must be right justified, and "F" formats must have the deci- 
mal point . The number of cards read per data block is listed 
at the bottom of each block. 

2 . Unformatted Data 

While the USER'S MANUAL data sheet defines SADX data 
in formatted fields of ten, the data may actually be read in 
a simplified fashion by separating data by commas or one or 
more blanks . If more than one number is contained on an un- 
formatted data card, a comma must appear somewhere on the card. 
If exponential numbers such as 2.5+10 are read on an unfor- 
matted card, there must be no imbedded blanks. Unformatted 
cards may be intermingled with formatted cards . Real numbers 
on an unformatted card must have a decimal point. 



33 



EXAMPLES 



Unformatted data 

5 . 7 .3 .2 .1. 3 + 6 ,-5 ,2 
Equivalent formatted data 

col 10 20 30 4^ 50 60 70 8£_ 

5 7 3.2 1.3+6 -502 

Unformatted data 

2 

2.3 
2 3 

Equivalent formatted data 

col 10 20 30 40 50 60 70 8^ 

2 

2 3 

2 3 

MOTE: The third line of data contains no commas and is 

therefore assumed to be already formatted. 

Placement of more than eight unformatted data on a 
card will create two (or more) formatted cards as required. 
Fields of zeros will be created if more data are required 
than are filled on an unformatted card. 

B. CONSTRAINTS 

Constraints are calculated and stored in the 0 vector 
as listed in the following chart. The total number of 
constraints 

NC0N = NFREQ + 2*MDSPLC+NLC*( 2*NEB + 4--NEF + NEB + 2*NEF) 
freq displ stress stress buckl buckl 
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Where NCON is the total number of constraints, NFREQ is the 
number of frequency constraints, NDSPLC is the number of 
displacement constraints, MLC is the number of load cases 
imposed, NEB is the number of bar or truss elements, and NEF 
is the number of frame elements . VJhen any of the constraints 
are missing from the G vector, all constraints are moved up 
in the vector. For example, if there is no frequency con- 
straint, then a displacement constraint would fill the first 
location of the G vector. 

CONSTRAINTS ARE STORED IN THE G VECTOR IN 
THE FOLLOWING ORDER; 

NFREQ FREQUENCY CONSTRAINTS 

2*NDSPLC JOINT DISPLACEMENT CONSTRAINTS 

STRESS CONSTRAINTS ARE STORED ELEMENT BY ELEMENT 

FOR A GIVEN ELEMENT CONSTRAINTS ARE STORED BY LOAD CASE 

FOR A GIVEN TRUSS ELEMENT AND LOAD CASE, 

CONSTRAINTS ARE STORED: 

TENSILE STRESS LOWER LIMIT 
TENSILE STRESS UPPER LIMIT 

EULER BUCKLING STRESS LIMIT (IF APPLICABLE) 

FOR A GIVEN FRAME ELEMENT AND LOAD CASE, 

CONSTRAINTS ARE STORED: 

NORMAL STRESS AT LOW NODE LOWER LIMIT 
NORMAL STRESS AT LOW NODE UPPER LIMIT 
NORMAL STRESS AT HIGH NODE LOWER LIMIT 
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NORMAL STRESS AT HIGH NODE UPPER LIMIT 

EULER BUCKLING STRESS LIMIT (IF APPLICABLE) 

LOCAL BUCKLING STRESS LIMIT (IF APPLICABLE) 

E . EXAMPLE 

The initial configuration of the braced cantilever beam 
is shown in Figure 5.1 Stress constraints are imposed as 
well as constraints on Euler and local buckling, displace- 
ment, and first fundamental frequency. A non-structural 
fixed mass is applied at the tip of the beam, and two load 
conditions (PI, and P2) are imposed. The structure's own 
weight will be considered as an imposed load as will be the 
fixed mass . 

The linking of design variables is demonstrated by link- 
ing the mean diameters of the two frame members (Dl, and 
D2 with D1=D2). Member size variables are; truss member 
areas (Al, and A2), frame member mean diameters (Dl, and 
D2 ) , and frame member thicknesses (Tl, and T2). Member 
sizing DESIGN VARIABLES are; XA(1)=A1, XA(2)=A2, XA(3)=D1=D2, 
XA(4)=T1, XA(5)=T2. 

Geometry variables are the attachment points (joints 4 
and 5) of the two truss members on the y and z axes (Yi^.ZS). 
Coordinate DESIGN VARIABLES are; XC(1)=Y4 and XC(2)=Z5. The 
structure's weight is then optimized with respect to member 
sizes and structure geometry for a total of seven design 
variables . 
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1 . Properties /Conditions 

Two material types are used: type 1, aluminum, is 

used for uhe truss members; and type 2, steel, is used for 
the frame members • The weight densities of the materials 
(p) are: 

type 1 = 0.1 Ib/in^ 

type 2 =0,3 Ib/in^ 

The Young's moduli of the materials (E) are: 
type 1 E = 10 .OE+6 psi 
type 2 E = 29. OE+6 psi 

The non-structural fixed mass attached at the tip of the 

structure weights 250 lb. The applied loads are: 

?1 1000.0 lb in the +y direction 

?2 1000.0 lb in the -z direction 

The acceptable maximum normal stresses are: 

tyre 1 -25000 psi < c < 25000 psi 

tyoe 2 -36000 psi < < 36000 psi 

^ — max — 

Displacement limits, imposed upon joint number 3 (the tip) 



each 


load 


case 


in the direction 


of loading. 


load 


case 


1 


y-direct ion 


+/- 3.0 in . 


load 


case 


2 


z-direct ion 


+/- 3.5 in. 



Bounds, placed on the positions of joints 4 and 5 along the 
y and z axes, are: 

joint 4 y-coordinate 0.0 inches to 200.0 inches 

joint 5 z-coordinate 0,0 inches to 100.0 inches 

The minimum natural frequency of the structure is constrained 
to be greater than IHz . 
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2 . Input Control Parameters 

The following input control parameters are given 



ease of 


following the 


example : 




NE3 = 2 


NEF = 2 


LO 

II 


NCJ = 3 


NMT=2 


IDVCLC=2 


NDJ = 2 


MEUBC=1 


LBUCK=1 


NFREQ=1 


NFMASS=1 


NEIG=1 


NEIG1=2 


NPRI=0 


NLC = 2 


MDSPLC=2 


NSTRES=2 


NSTW=1 


NFMW=1 





Table VI is a listing of commonly used nomenclature. 

The following USER'S MANUAL is divided into blocks 
A through P. Appearing directly below each data field line 
are the parameters for the TRUSS-BRACED CANTILEVER BEAJi 
example . It is important to note that the user may choose 
any consistent system of units . 



38 



TABLE VI 



COMMON VARIABLE NOMENCLATURE 



A 

BL 

BO 

CHASDIM 1, 2- 
DIH 



EPSEIG 

FC 1 , , . . FC N- 

FX 

FY 

FZ 

GR A? 

ID VC LC 



IX 

lY 

IZ 

IX X 

lYY 

IZZ 

JN 

LBOCK 

LC 

L.1ASS 



LSO 
L3 SCT 

SCJ 
ND J 

ND SG 1 
NDSG2 
NDSG3 

NDSGU 



MEMBER’S CROSS-SECTIONAL AREA 
LOWER BOUND ON DISPLACEMENTS 
UPPER BOUND ON DISPLACEMENTS 
CHARACTERISTIC DIMENSIONS FOR FRAME MEMBERS 

FOR LSECT.EQ.1: MEAN DIAMETER AND THICKNESS 
DIRECTION 1=X, 2=Y, 3=Z, 4=rot about x, 

5=rot about y, 6=rot about z. 

YOUNGS MODULUS 

CONVERGENCE TOLERANCE OF EIGENVALUE 
SOLUTIONS f DE?AOLT=.0001) 

LOWER BOUND ON FIRST, SECOND , ETC. 

SYSTEM NATURAL FREQUENCY 
LOAD FORCES APPLIED IN THE X DIRECTION 
LOAD FORCES APPLIED IN THE Y DIRECTION 
LOAD FORCES APPLIED IN THE Z DIRECTION 
ACCELERATION DDE TO GRAVITY fDEFAULI= 386. U) 
THE OPTIMIZATION OPERATION IDENTIFIER 

1 FOR OPTIMUM MEMBER SIZE ONLY 

2 FOR BOTH OPTIMUM MEMBER SIZE AND GEOMETRY 

3 FOR OPTIMUM GEOMETRY ONLY 
CONSTRAINT IDENTIFIER. IF NON-ZERO THE X-DO? 

IS CONSTRAINED 

CONSTRAINT IDENTIFIER. IF NON-ZERO THE Y-DCF 
IS CONSTRAINED 

CONSTRAINT IDENTIFIER. IF NON-ZERO THE Z-DCF 
IS CONSTRAINED 

CONSTRAINT IDENTIFIER. IF NON-ZERO THE 

ROTATION ABOUT THE X-AXIS IS CONSTRAINED 
CONSTRAINT IDENTIFIER. IF NON-ZERO THE 

ROTATION ABOUT THE Y-AXIS IS CONSTRAINED 
CONSTRAINT IDENTIFIER. IF NON-ZERO THE 

ROTATION ABOUT THE Z-AXIS IS CONSTRAINED 
JOINT NUMBER (GLOBAL) 

KEULEP- EULER BUCKLING COEFFICIENT FOR BAR 
ELE ME NTS 

LOCAL BUCKLING CONSTRAINT IDENTIFIER 

IF LBUCK.NE.O, LOCAL BUCKLING CONSTRAINTS 
WILL BE APPLIED TO TUBULAR FRAME MEMBERS 



LOAD CONDITION 
LUMPED MASS OP 
IF LMASS.EQ 
IF LMASS.NE 
ELEMENT NUMBER 
CROSS-SECTION 
LSECT.EQ. 1 
THE NUMBER OF 
THE NUMBER OF 
DESIGN VAR I 
TRUSS MEMBER A 
COORDINATE DES 
FRAME MEMBER C 
DESIGN VARI 
FRAME MEMBER C 
DESIGN VARI 



TIONS 

.0 CONSISTENT MASS MATRIX USED 
.0 LUMPED MASS MATRIX USED 



TYPE IDENTIFIER 
INDICATES A TUBULAR MEMBER 
JOINTS 

(OR LINKED TO) 



CONSTRAINED 
JOINTS WITH 
ABLE S 

HSa“ DESIGN VARIABLE NUMBER 
IGN VARIABLE NUMBER 
HAEACTBRI3TIC DIMENSION 1 
ABLE NUMBER 

HARACTSRI3TIC DIMENSION 2 
ABLE NUMBER 
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NDSPLC 

NEB 

NEF 

NEIG 

NEIG1 

NEUBC 



NFHASS 

HFMW 



NFHEQ 

NID 

NJ 

NLC 

NLJ 

NMT 

NPRI 



NST5ES 



NSTW 



FOISSN 

FHO 

SIGMIM 

SIGMAX 

TX 

TY 

TZ 

xa 

XAL 
XA U 
XC 
XCL 
XC'J 



- THE NUMBER OF DISPLACEMENT CONSTRAINT SETS 

- THE NUMBER OF TRUSS ELEMENTS 

- THE NUMBER OF FRAME ELEMENTS 

- NUMBER OF PRECISE EIGENVALUES TO BE FOUND 

- NUMBER OF EIGENVALUE/EIGENVECTOR TD BE 

EVALUATED DEFAULT =MIN{2*NEIG, NSIG+8) 

- EULER BUCKLING CONSTRAINT IDENTIFIERS. 

IF NEUBC. NE.O EULER BUCKLING CONSTRAINTS 
WILL BE IMPOSED ON BAR AND FRAME ELEMENTS 

- NUMBER OF NON- STRUCTURAL FIXED MASSES 

ATTACHED TO THE STRUCTURE 

- FIXED MASS WEIGHT IDENTIFIER 

IF {NFMW. NE. 01 FIXED MASSES MILL BE 
CONSIDERED AS LOADS 

- NUMBER OF FREQUENCY CONSTRAINTS 

- NUMBER OF INDEPENDENT DEGREES OF FREEDOM 

- THE NUMBER OF JOINTS 

- NUMBER OF LOADING CONDTIONS 

- NUMBER OF LOADED JOINTS FOR THIS LOAD CONDITION 

- NUMBER OF SEPARATE MATERIAL TYPES 

- INPUT PRINT CONTROL 

IF NPRI. NE.O INPUT VALUES 3iT~ NOT PRINTED 
IF NPRI.EQ.5 LOCI/LOCR/IA/P.A ARRAYS PRINTED 

- THE POPCE/STRESS PRINT IDENTIFIER 

IF NSTRES. EQ.O MEMBER FORCES PRINTED 
IF NSTRES. EQ.1 MEMBER STRESSES PRINTED 
IF NSTRES. EQ. 2 BOTH are PRINTED 

- STRUCTURE WEIGHT IDENTIFIER 

iF (NSTW. NE. 0) THE STRUCTURE’S OWN HEIGHT 
HILL BE CONSIDERED AS LOADS 

- POISSON’ S RATIO 

- MATERIAL WEIGHT DENSITY 

- MINIMUM ALLOWABLE NORMAL STRESS 

- MAXIMUM ALLOWABLE NORMAL STRESS 

- TORSIONAL MOMENT APPLIED ABOUT THE X-AXIS 

- TORSIONAL MOMENT APPLIED ABOUT THE Y-AXIS 

- TORSIONAL MOMENT APPLIED ABOUT THE Z-AXIS 

- INITIAL VALUE OF AREA DESIGN VARIABLE 

- LOWER BOUNDS 0 N XA 

- UPPER BOUNDS ON XA 

- INITIAL VALUE OF COORDINATE DESIGN VARIABLE 

- LOWER BOUNDS D N XC 

- UPPER BOUNDS D N XC 
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DATA BLOCK A 



DESCRIPTION: Title Card 



Format ar.d Example 



TITLE 



FOR HAT 



20A4 



TRUSS-BRACED CANTILEVER BEAM 



F IEL D CO NTEN TS 

1 ANY 80 CHARACTER TITLE MAT BE GIVEN ON THIS 
LINE 



4 





■X 



I 



I 



MIh. lliOCK 5 

DE SCS I? TIP K: Control Parameters 
Format and Example 



NEB 


NEF 


NJ 


NCJ 


NMT 


IDVCLC 


NDJ 


FORMAT 
















7110 



2 


2 


5 


3 2 2 


2 














NEUEC 


LBUCK 


NFREQ 


NFMASS 


NEIG 


NETG1 


NPRI 


f ormar 
7110 












1 


1 


1 


1 


1 


2 


0 



NLC 


NDSP LC 


NS TEES 


NSTN 






FORMAT 

5110 












2 


2 


2 


1 


1 



MOTE: DEFINITIONS OF PROGRAil INPUT CONTROL 
PARAMETERS ARE LISTED ON NEXT PAGE 
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I 

s 



I 



^ I 

t 



t 



1 1 



F IEL D C ONT ENT 

FIRST CARD 

1 NEB-number of bar alements 

2 NE?-number of frame sleraenrs 

3 NJ-number of joints 

4 NCJ-number of constrained joints 

5 NHT-number of seperate material types 

6 IDVCLC-design variable control parameter 

If (IDVCLC. EQ. 1) NDV=NDVAR1 
optimizes wrt member size only 
If (IDVCLC. EQ. 2) NDV=NDVAR1+ NDVAR2 
optimizes wrt member size and geometry 
If (IDVCLC. EQ. 3) NDV=NDVA82 
optimizes wrt geometry only 

7 NDJ-total linked and design variable joints 

(i. e. number cf ’movable’ joints) 

SECOND CARD 

1 NEUBC-Euler buckling constraint identifier 

If (NZUBC.NS.O) EULER buckling constraints 
will be imposed on bar elements 

2 LBUCK-local buckling constraint identifier 

If (LBUCK.se. 0) local buckling constraints 
will be imposed on tubular frame elements 

3 NFREQ-number of frequency constraints 

4 NFM ASS-nu mber of fixed masses attached to 

structure 

5 NEIG-number of precise eigenvalues to 

be evaluated 
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6 NEIGI-number of eigenvalues to be evaluated 

DS?AOLT=min. of (2*NSIG , NEIG+3) 

7 NPRI-prinx control identifier 

NPRI.ne.O inpux info not printed 
NPR1.eq.5 RA/ lAAOCR/LDCI arrays 
will be printed for debugging 

THIRD CARD 



1 


NLC-n 


umber of 


1 oad 


cond 


it : 


Lo 


ns 




2 


NDSPL 


C-nu raber 


of 


displa 


came 


nt consrr 


ain 


3 


NSTRE 


S-fcrcs/ 


St r 


ess pri 


nt 


1 


dent if ier 






Tf 


{NSTRES. 


EQ . 


0) 


stre 


3 3 : 


zS 


are prin 


t 9d 




If 


(N ST EES. 


SQ . 


1) 


fore 


3 5 


a 


T9 t)rir.t3 


d 




If 


(NSTRES. 


3Q . 


2) 


both 


1 


- 


print 9d 




4 


KSTW- 


structur 


e w 


eight i 


de: 


It 


ifier 






If 


(NSTW.N3 


.0) 


rhe St 


ruct 


ure's wai 


ght 



considered as loads 

5 NFH»-fixed mass weight iianrifier 

If (HFHW.MS.O) the fixed masses are 
considered as loads 
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DATA 5L0« C 

DESCRIP TION; Dynamic Analysis Information 
Format and Example 



LflASS 


GRAY 


EPSEIG 




FOSHAT 
I10,2?10. 0 


0 


386.4 


o 

• 

o 






FIELD 











1 LMASS- lumped mass option (if LSASS.N3.0) 

the lumped mass matrix is used. 

2 GRAV-accleration due to gravity 

(defa ult=3 86 . 4 inches/sec^) 

3 SPSEIG-convergence tolerance on eigenvalue 

solution, (de faul"=.000 1) 





I 



I 



I 



DATA BLOCK D 



DE SCR I PTION : Joint Coordinates 

Forma- and Example 




FIELD CON TENTS 

1 JN- joint cooriinane number 

2 X-x coord inane 

3 Y-y coord inane 

a z-z coord inane 



NOTE: Number of cards raai=NJ 
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Mil block s 

Omit this block, if SDJ=0 in BLOCK B 
DE SCR IP TIP N: Coordinate Design Variable Linking Data 

Format and Example 



JN 


IX 


lY 


IZ 


PCX 


PCI 


PCZ 


FO.RMAT 
4110, 3f 10,0 


















4 


0 


1 


0 


1.0 


1.0 


1 . 0 
















5 


0 


0 


2 


1, 0 


1.0 


1 . 0 



FIELD 



CONTENTS 



1 


IX- design variable 


associ ared 


with 


X 


coord 


2 


lY-design variable 


associat e,d 


with 


7 


coord 


3 


IZ-design variable 


associated 


with 


Z 


coord 


4 


PCX-participati on 


coefficient 


of X- 


c 


oord . 



5 ?CY-narticipati cn coefficient of y-coord. 

6 PCZ-oarticioati on coefficient cf z-coord. 



raad=NDJ 



NOTE: Number of cards 





I 




D ATA BLOCK F 



DE SCB IP TIO N: Material Properties 
Format and Example 



E 


SHO 


SIGHIN 


SI GHAX 


KEULER 


POISSN 


FORMAT 
6F10. 0 
















1. OE+7 


0. 1 


-250 00, 


25000. 


2.0 


0.27 














2.9E+7 


0.3 


-360 00. 


3 60 00 . 


2.0 


0.27 



FIE LD COM TE NTS 

1 E-Young*s Modulus 

2 3H0-marerial density 

3 SIGMIN-mi nimum allowable stress 

4 siGNAX-maximum allowable stress 

5 KEULSR-Euler buckling coefficient 

6 POISSN-Poisson* s ratio 



MOTE: Number of cards rsad=NHT 
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I 











NS3=0 in 3L0CK 3 



SilA block G 

Omit this block if 
D 3SC 3IP TI0N : Bar Element Information 
Format and Example 



LNC 


NODE1 


NODE2 


MAT COD 


NDSG1 


A 


LSECT 


FORMAT 
5110 
f 10,110 


















5 


2 


4 


1 


1 


3.0 


1 
















2 


2 


5 


1 


2 


3. 0 


1 



FIELD C ONT ENTS 

1 LMO-elemenr number 

2 NODEl-glcbai number associated wich the 

elemenr’s 1st node 

3 NODE2-global number associaned with the 

element's 2nd node 

U h ATCOD-nater ial type of this element 

5 NDSGI-desian variable number associared with 
this element's area 

6 A-member cross-sectional area 

7 LSECT-cro ss-ssction type identifier 

LSECT.EQ.1 indicates a tubular member 

NOTE: Number of cards read=NEB 
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DATA BLOCK H 

Omit this block if N3F=0 ir. BLOCK 5 
D ES CRI PTION ; Frame Element Information Forma* and 
Example 



LNC 


NODE1 


NODE2 


MAT COD 


NDSG3 


NDSG4 


LSECT 


FORMAT 

7110 
















3 


1 


2 


2 


2 


4 


1 






CHAHDIM1 


CHARDIM2 




FORMAT 

2F10 



5.0 


1.0 




u 


2 


3 


2 


3 


5 


1 


o 

f 

in 




1.0 





FIEL D S2EIM2S 



1 LNO-element number 

2 NODEI-global number associated winh the 

element's 1st node 

3 N0DS2-global number associated with the 

element's 2nd node 

4 MATCOD-ma ter ial type of this element 

5 NDSG3-design variable number associated with 

the element's 1st characteristic dimension 

6 NDSG4-design variable number associated with 

the element's 2nd characteristic dimension 

7 LSSCT-crcss-section type identifier 

LSECT.SQ.1 indicates a tubular member 
NOTE: Number of cards read=NSP 
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A 



I 



4 



DA^A BLOCK I 



D 5SC HI PTICN ; Joint Constraint Information 
Fermat and Fxampls 



JN 


IX 


IT 


IZ 


I XX 


lYY 


IZZ 


FOaXAT 
















7110 



1 


1 


1 


1 


1 


1 


1 




2 


0 


0 


0 


0 


0 


0 




3 


0 


0 


0 


0 


0 


0 




4 


1 


1 


1 


1 


1 


1 




5 


1 


1 


1 


1 


1 


1 



■p T-S> T n 

■ 4 . ^ ^ XJ 



CONTENTS 



1 

2 

3 

4 

5 

6 
7 



JN- joint number 

IX- x-dis place ment constraint, iaentii 
IT- y-dis placament constraint idsntif 
IZ- z-displacenent constraint identif 



IXX- x-axis 


rot 


arion 


constraint 


i d an 


ITT- y-axis 


rot 


atior. 


constraint 


idea 


IZZ- t-axis 


rot 


atior 


constraint 


iden 



if . NS. 0 - CO rresponding DO? constr 



f ier . 
fier . 
ined 



NOTE: Number of cards raad=NCJ 
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slock j 

Omit this block if flLC=0 in BLOCK 3. 
DE S CR I?T IO N: Joint Loading Information 

Format and Example 



NLJ 




FORMA? 






no 



1 



JN 


FX 


FY 


F2 


TX 


TY 


TZ 


FORMAT 
1 10,6F10 




















2 


o 

• 

o 


1000. 


0.0 


O 

« 

o 


o 

• 

o 


o 

• 

o 




1 


















o 

• 

o 


o 

• 

o 


-1000 . 


o 

• 

o 


o 

• 

o 


0.0 





FIELD COHTENI 

1 NLJ-number of loaded joinrs for this load 

condition 

1 JN- joint number 

2 FX 

3 FY- Forces in the X,Y,Z directions 

4 FS 

5 TX 

6 TI- aomenrs about the X,Y,Z axes 

7 TZ 

NOTE: Number of cards read per sei:=NLJ 
Number of sets of cards rsad=NLC 
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« 



I 






DATA BLOC K K 

Omit this bloclc if NFMASS=0 ia BLOCK B 
DE SCR IPTION: Fixed Mass Information 

Format and Example 



JN 


MASS 




FORMAT 

no, F 10 








3 


250. 0 



FI3L D C ONT ENTS 

1 JN- joint number 

2 MASS-poinr mass at joint (JN) in force unirs 



NOTE: Number of cards read=NFMASS 




I 



I 



DATA BLOCK L 

Omit this block if IDVCLC=3 

DESCRIPTION; D'ssign Variabls Informa tioTi 

(MEMBER SIZE VARIABLES) 

Format and Example 



XA (1) 


XA (2) 


.... XA(ND7AR1) 


FORMAT 

8F10.0 










20. 0 


20.0 


20.0 


2.0 


2.0 












XAL (1) 


XAL (2) 


. . . .XAL(NDVARI) 


FORMAT 

8F10.0 










0.50 


0.50 


4.0 


0.10 


0.10 












XAU (1 ) 


XAU (2) 


. . . .XAU(NDVARI) 


FORMAT 
8F1 0. 0 










25.0 


35.0 


2 5.0 


2.5 


4.0 





FIELD CONTEN TS 



1 XA- initial value of area design variables 

2 XAL-lower bounds on area design variables 

3 XAU-upper bounds on area design variables 

MOTS; read one value of XA, XAL, XAU for each 
independent member size variable defined in 
BLOCKS G and H 

Number of cards read = as required 
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DATA BLOCK M 

Omit this block if IDVCLC=1 

DES CRIP TION; Design Variable Information 

(COORDINATE VARIABLES) 

Format and Example 



XC (1) 


XC (2) 


.... XC(NDVAR2) 


FORMAT 








3F1 0. 0 



150.0 


5 0.0 








XCL (1 ) 


XCL (2) 


. . . .XCL(NDVAR2) 


FORMAT 

3F10.0 








0.0 


O 

• 

o 






XCO (1) 


XCU (2) 


. . . .XC0(NDVAR2) 


FORMAT 

3F10.0 








2 00.0 


100.0 



FIELD CONTENTS 



1 


XC- initial value 


cf 


coord. 


design variaoles 




XCL-lower bounds 


on 


coord. 


design variables 


3 


XCO-upper bounds 


on 


coord . 


design variable? 




NOTE: read one va 1 


ue 


of XC,X 


CL, XCU for each 



independent coord, variable defined in BLOCK D 
Number of cards read =as required. 
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a 




J 



i 





D ATA BLOCK N 

Omir this block if NDSPLC=0 in BLOCK B 
DSSCRI PTION: Joint Displacement Constraint In formation 
Format and Example 







CON TENTS 



1 JN- joint number 

2 DIS-direction 

1=x,2=y,3=z displacement 
4=x,5=y»6=2 rotation 

3 LC-load condition 

4 3L-lower bound on displacement 

5 3U-upper bound on displacement 



NOTE: Number of cards read= NDSPLC 
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I <1 

I 




D ATA BLOC K 0 

Omit this block if MFRBQ=0 in BLOCK 3 
DESCR IP TIO N; Frequency Cor.sc rein- In f orma-icr. 
Format and Example 



?C1 


FC2 


FC3 




FORMAT 










8F10.0 



1.0 



field con 21 NTS 



1 


?C1- lower bound on 


first natural frequency 




constraint, in Hz, 


(cycles per second) 


N 


FCN- lower bound on 


NF3E2-th natural 




frequency ccnstra 


int in Hz. 




(cycles per 


second) 



NOTE: Number of cards read = as required 



MIk BLogs ? 



D3 SCE IP TI0 N: Snd card 
Format and Example 




NOTE: This card ausT appear ar rhe end of the 



SADX data. 
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VI . NUMERICAL EXAMPLES 



A. INTRODUCTION 

Design of three-dimensional truss and frame structures 
are presented herein and the corresponding numerical results 
are summarized to demonstrate the use of the SADX code. 

The examples begin with the truss-braced cantiliver 
beam. 

B. EXAMPLE 1: TRUSS-BRACED CANTILIVER BEAM 

The simple truss-braced cantiliver beam, as shown in 
Figure 5.1, has been previously used for the SADX USER'S 
MA.NUAL example. It consists of two steel tubular frame 
members with a common outer diameter and different wall 
thicknesses arranged as a cantilever beam along the X-axis . 
There is a fixed 250 pound mass at the tip of the beam. 

Two aluminum truss members are attached from the beam .mid- 
point to points on the Y and Z axes . This structure is 
designed for optimum me.mber size and geometry under a set 
of two load conditions and subject to constraints on maximum 
stress, maximum joint displacement, Euler and local buckling, 
and minimum structure natural frequencies. The weight of 
the non-structural fixed mass and the structure's own weight 
are i.mposed as loads. There are five member size and two 
coordinate design variables, and a total of 41 constraints. 
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Figure 5.1 TRUSS -SR^.CED CANTILEVER BEAM 
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ihe number of analyses required for this design is 116 using 
2.75 seconds of CPU time and terminating on the 12th itera- 
tion. Or the analyses conducted, 84 were recuired for the 
calculation or gradients. The weight of the structure, 
including non-structural fixed masses, is reduced from 9542 
pounds to 333 pounds. Results are given in Table VII. 

C. EXAMPLE 2: T.i70-TIER 3-D PORTAL ?R.^.ME WITH TRUSS X-3RACES 

The two-tier three-dimensional portal frame with truss 
member diagonal braces, as shown in Figure 6.2, is a symme- 
tric moderately sued case that can be analyzed easily by a 
variety of other codes. There are 20 truss elements, 16 
frame elements , and four non-structural fixed masses . The 
material used is steel. This structure is designed for 
opti.mum member size and geometry under a set of three load 
conditions and subject to constraints on maximum stress , 
maxim’im joint displacement, Euler and local buckling, and 
minimum structure natural frequencies . The weight of mhe 
non-strucrural fixed masses and the structure's own weight 
are imposed as loads . There are 10 member size and five 
coordinate design variables and a total of 493 constraints. 
The number of analyses required for this design is 335 using 
120 seconds of CPU time and terminating on the 21st itera- 
tion. Cf the analyses conducted, 305 were required for the 
calculation of gradients. The weight of the structure, 
including non-structural fixed masses, is reduced from 9302 
pounds to 1462 pounds. Results are given in Table VIII. 



51 



PIN4L CPTIKIZATICN INFORMATION 



Cej = C.5675CTS+C3 

TKEFE ARE 1 ACTIVE CONSTRAINTS 

CChETRAINT NUMBEFS APE 
36 

ThEFE ARE 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

TESMINA'^ICN CRITEPICN 

AES (C3J ( I I-G3J ( I-l n LESS THAN CABFUN FCR 3 I TE 
NLMBEP CF ITcRATICNS » 12 

C5JECTIVE FLINCTICN WAS EVALUATED llA 

CCNSTRAIN7 FUNCTIONS WERE EVALUATED llA 

THIS RUN FECUIREC 116 STRUCTURAL ANALYSES 

NUMEER OF SECCNCS REQUIRED FOR EXECUTION = 2.79 

WEIGHT OF STRLCTUFE GIVEN AREAS S LENGTHS 
WiIGHT= 0.5S751E+C3 



TOTAL WEIGHT INCLLDING FIXED MASSES 
TC7AL W£ICHT= C.S3751E+03 



JOINT COCFDINATSS 
JCINT X 

2 cIl(CCOH+C3 

3 O.2C0Ca5+C3 

4 0,0 

5 0.0 



0.0 
G. 0 
0.0 

C. 54259E+02 

C. 0 



0 .0 
0.0 
0 .0 
0 .0 

0.^5432E+02 



ELEMENT 

1 

2 



ELENENT INFCRMATICN FCri 2^R ELEMENTS 
cLEMENT-JCIM F cL JiTI CNSHI ?S 
1 NOCE 2 
5 

ELE^^ENT INFCRvaTICN FCR FRAME ELEMENTS 
ELE^ENT-JCINT R EL^a TI CNSHI PS 



NODE 

2 

2 



AREA 

0.3^8 lE*f-C 1 
0. 3499E^01 



length 
0. 13C3E>03 
0. 1098 E +03 



LNC NCCEl NCCE2 LENGTH 

3 1 E 0.3329E+01 O.IOOOc+CB 

423 0.236£E+C1 0.10COE+C3 



CHAR.OIM.l 
0 .Ac39E+01 
0.4329E+C1 



EIGENVALLES ANC E IGcNYECTCRS 




VCCE 


NL'HEER 1 






FRECU 


ENCY = C.^1=8EE+C2 CPS 






C .3631 


9E+C2 RPS 




eigenvalue = C.6927CE+05 




EIGENVECTOR DEGREE 


CF FREEDOM 




JCINT 


1 


2 


3 


1 


c.o 


0.0 


0*0 


2 


c.o 


c. c 


0. 0 




•0 .224253^01 


-C.32985E-01 


-0. 745C76-01 


4 


-0.50777E-C1 


-C. 107^5 E-01 


0.15330E-03 


c 


-0.23435E-C1 


-C. 1T155E-01 


C. lOOCCE+01 



RATIONS 

TIMES 

TIMES 



CHAR. DIM. 2 
0.2190E+00 
0. 155EE+00 



example 2 



Y 




Figure S.2 TWO-TIER 3-D PORTAL FRAME WITH TRUSS X-3RACES 
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EXAMPLE 2 : 



EIMAL OPTIMIZATION INECRMATIOM 



FIN/L CFTiyiZATICN INFORMATION 

C8J » C.1^4ie3E+C4 

TFEF5 A<?£ 3 ACTIVE CONSTRAINTS 

CONSTRAINT NUMBERS ARE 
A6 SI ISA 



THERE ARE 



C WCLATED CCNSTRAINTS 



THERE ARE 1 ACTIVE SIDE CONSTRAINTS 
CECISION \APIAeLES AT LCWSR CR UPPER BOUNDS 
(MINUS INCICATES LOWER SOUND) 

-14 

TERNINATICN CRITERION 

AES (C5J ( I 1-CBJ ( I-l n LESS THAN CAB RJN FCR 3 ITERATIONS 
NUMBER OF ITERATIONS = 21 

OBJECTIVE FLNCTICN WAS cVALUATEC 333 TIMES 

CONSTRAINT FUNCTIONS WERE EVALUATED 3B3 TIMES 



THIS RUN PECUIREC 385 STRUCTURAL ANALYSES 

NUMBER OF ScCCNCS RECUIREC FOP EXECUT ICN » 120.7 

WEIGHT CF STRUCTURE GIVEN AREAS S LENGTHS 
WEIGHT* 0.1441EE+C4 



TOTAL WEIGHT INCLL'CING FIXED MASSES 
TCIAL WEIGHT* C, 16418E+04 



JOINT COCRDINATES 



JCINT 


X 




Y 


Z 


1 


0.6 3385 


E-^C2 


0.0 


0 .9C0C2E+02 


2 


0 . AE2e5£+ C2 


G. 0 


-0 .90002E-i^02 




-0. 8 -385 


E^C2 


C.O 


-0.9C0C2E+02 


4 


-0.63335 


E*C2 


0.0 


0.9C0C2E-»^02 


5 


0 . aC239E-»-C2 


C.20C00E*02 


0 .69074E+02 


6 


0,aC239 


E+C2 


C.20000E+C2 


-0.69074E+02 


7 


-0. 8C239 E+C2 


C.20000E+02 


-0 .69074E+02 


a 


-0. 6C239 


E+C2 


C. 20000E+02 


J.89074E-I-02 


9 


0. 1C0C0E + C3 


C. 10000E-KC3 


0. 1C0CC5+03 


1C 


O.ICCCO 


E-^C3 


O.lOOOOE+03 


-0 .lCOOOE+03 


1 1 


-0. 1CCC0E+C2 


C. 100C05 + 03 


-o.icocoe+03 


12 


-0. ICCCO E+C3 


C.100006+C3 


0 .lCOCOE+03 


E L E V E NT 


INFCRXATICN FCR 


3AR ELEMENTS 






JCINT PEL 


ICNS 


pIPS 




ELEMSNT 


NODE 1 


NCCE 2 


AREA 


LENGTH 


1 


1 


6 


0.^255E+C0 


0. 16125+03 


2 


2 


5 


0.4255E-^CO 


0. 1612E+03 


3 


2 


7 


0. ^2 5 5E + Cu 


0. 1465E + 03 


4 


*2 


6 


0.4255E+C0 


0. 1465E+03 


5 


2 


A 


C.^255c-t-C0 


0. 1612E+03 


6 


4 


7 


0. 4255E+C0 


C. 1612E + 03 


7 


4 


5 


0. *+2 55E-^C0 


0. 14c5£+03 


a 


1 


a 


C. 4Z55E+C0 


0. 1465E+03 


9 


e 


7 


0. 4255E+C0 


C • 21135+03 


10 


6 


3 


0 • 4255 E-^CO 


0. 2113E+03 


11 


5 


IC 


0. 1544E-^Ci 


0. ISSlE+03 


12 


6 


9 


C. 154-i«f= + 01 


0. 1381E+03 


*2 


6 


11 


0. 1544£^^C1 


0. 1996E+03 



D. EXAMPLE 3: DD-963 FOREMAST 



Example three is a redesign of the forward mast on the 
DD-963 of SPRUANCE class destroyer. The DD-963 foremast has 
been chosen as the third test case for the following reasons : 
1) the structure is typical of the masts found on many com- 
batants in the United States and other navies, 2) high top- 
side weight reduction is desirable from a stability viewpoint 
for any ship, 3) the structural members are predominantly 
tubular, 4) the problem is sufficiently complex to make 
conventional design methods cumbersome, 5) the member size 
and loading information is available. 

The structure as shown in Figures 6.3 through 6.5 is 
constructed of 172 aluminum frame members. The material 
used in the structure is 5086 H32 aluminum, an alloy with 
moderate strength, good weldability, and good corrosion 
resistance . 

This structure is designed for optimum member size under 
a single load condition and subject to constraints on maxi- 
mum me.mber stress, maximum joint displacement, Euler and 
local buckling. 

Some structural simplifications are made. The weights 
of mast-mounted radars, antennas, and platforms are modeled 
by a series of fixed masses which are imposed as loads along 
with the structure's own weight. The forces due to halyards 
and wire antennas are applied as loads . Inertial forces due 
to ships motion are calculated for the initial design point 
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Figure 6.3 DD-363 FOREMAST PORT SIDE 
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5 7 
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and consisu of The prcducT of a nerier’^s half-weiThT and a 
coondinaT 6 -dependenf load-facTor aoolied an each end of The 
elemenT [Ref. cl. The load-facTor milTiolier is deoendenT 
on The dis Tances f roni a ooinT aT The inTersecTion of “he 




INC 

1 

2 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

1 ^ 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

3C 

31 

32 



TABLE IX 



EXAMPLE 3: FINAL OPTIMIZATION INFORf^A.TIOM 



FINAL OPlIf-I ZA7ICN INFCRMATICN 

CEJ = C.3986C4E-^05 

2 ACTIVE CONSTRAINTS 
CCNSTRAINT NUMEERS ARE 
660 666 

THERE ARE C VICLATEO CONSTRAINTS 

THERE ARE C ACTIVE SIDE CONSTRAINTS 

TEF^^INATICN CRITEPICN 

ABS(08j( n-o0j(i-i) less than cabrjn for 

NU^'BER CP ITERATIONS = 17 

OBJECTIVE FUNCTION WAS EVALUATED 
CCNSTRAINT FUNCTIONS WERE EVALUATED 



3 ITERATIONS 

553 TIMES 
553 TIMES 



THIS RUN FEOUIREC 555 STRUCTURAL ANALYSES 

WEIGHT OF STRUCTURE GIVEN AREAS L LENGTHS 
WEIGHT^ 0.29860E+05 

TOTAL WEIGHT INCLUDING FIXED MASSES 
TOTAL WEIGHT^ 0.56746£-f05 

NUMBER OF SECONDS RECUIRED FOR EXECUTION = 4482.46 

NC. OF FFAME ELEMENTS ^ 172 

ELEMENT-JOINT RELATIONSHIPS 



NCCEl 


N0CE2 


AREA 


LENGTH 


CHAR .DIM.l 


CHAR .DIM. 2 


4 


12 


G.5979E+02 


0. 1061 E+C3 


C.19CC5+02 


0.1002E+01 


-a 


11 


0. 5979E+C2 


0. 1061E+03 


C.19C0E+O2 


0.1CC2E+01 


1 


c 


C. 59T9£-t-C2 


0.1112E+03 


C. 19005+02 


0.1002E+01 


2 


10 


C.59T9E>C2 


0. ill2E+03 


C. 1900E+02 


0.10C2E+01 


4 


8 


0.62475^01 


0.1342E+03 


C. 65635+01 


0.203CE+00 


3 


6 


C.62^7E+C1 


0. 1342E+03 


C.6563E+01 


0.30306+00 


1 


a 


0.62^7E+01 


0. 1369E+02 


C. 6563E+01 


C.2C3CE+00 


2 


6 


0.6247E+01 


0.13696+03 


0.6563E+01 


0.3030E+00 


8 


12 


C.62^7E>C1 


0. 1295E + 03 


C.6563E+01 


0.3030E+00 


6 


11 


0.62A7E-*'01 


0.1295E+03 


C.O563E+01 


0.303CE+00 


6 


c 


C.62A7E+01 


0.1132 6+03 


C.6563E+01 


0.3030E+00 


6 


1C 


0.62A7E+C1 


C. 1132E+03 


C.6563E+01 


C.3C30E+00 


3 


7 


0.6247E-+01 


0.6846 E+02 


0.65635+01 


0.303CE+00 


4 


7 


C.6247E+C1 


0.6846E+02 


C.6563E+01 


0.3020E+00 


7 


11 


0.6247E^C1 


0.6846E+02 


C.6563E+C1 


C.3G3CE+00 


7 


12 


0.62^75+01 


0.6846 E+02 


0.6563E+01 


0.3030E+00 


1 


c 


G.6247E+C1 


0. 1196E+03 


0.6563E+O1 


0.3030E+00 


2 


c 


0.62^7E+C1 


0.1196 5+03 


C.6563E+01 


0.2G3CE+00 


5 


9 


0.62475+Cl 


0.1097 E+03 


0.6563E+01 


0.3030E+00 


5 


1C 


C.62^7E+C1 


0. 1097E+03 


C. 65635+01 


0.3C3CE+00 


9 


22 


0.6247E+01 


0.21726+03 


0 .6565 5+01 


0 .2030E+00 


1C 


11 


C. 6247E+C1 


0. 21726+03 


C. 65635+01 


0.3020E+00 


11 


12 


0.6247E+01 


0.84COE+02 


C. 6563E+01 


C.3C3C6+00 


c 


1C 


0.62475+Cl 


0.1894E+03 


0.65655+01 


0.30306+00 


12 


21 


0. 5979E+C2 


C. 1061E+03 


C. 19C0E+O2 


0. 1002E+01 


11 


16 


0. 5979E+02 


0.1081 E+03 


G.1900E+02 


0. I0C25+01 


9 


16 


C.5979E+C2 


0.1111 E+03 


C. 19005+02 


C.1002E+01 


10 


17 


C.5979E+G2 


O.llllE+03 


C. l<500E+02 


0. 1CC2E+01 


9 


21 


0.7358E+C1 


0.2380E+03 


0.77o7E+01 


0.3C16E+00 


1C 


16 


C.7356E+C1 


0.2380E+03 


C.7767E+01 


0.3016E+00 


11 


12 


0.7356E+01 


0. 6846E+02 


C. 7767E+01 


0.3C16E+00 


12 


13 


0.7358E+01 


0.68466+02 


C.7767E+01 


0 . 30 16E+00 
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1 








I 






13 

13 

1C 

u 

17 

18 

16 

12 

11 

20 

18 

IS 

21 

18 

16 

17 

16 

17 

20 

1 S 

20 

IS 

18 

21 

2A 

24 

16 

17 

26 

27 

29 

28 

31 

30 

28 

15 

15 

30 

29 

28 

26 

27 

26 

27 

29 

28 

26 

27 

28 

29 

33 

* 2*3 

34 

36 

39 

37 

39 

37 

34 

36 

34 

36 

40 

38 

25 

25 

26 

27 

34 

3 5 

35 

35 



18 


0.7358E+C1 


0.6846 5+02 


0. 77675+01 


2 1 


c.735ee<-ci 


0. 6846E+02 


C. 77675 + 01 


16 


0.7358E+01 


0.2094E+03 


0.77675+01 


21 


C. fA2^e + Cl 


0.1866 5+03 


0.5024E+01 


18 


0. 5^2AE+C1 


C. 1866E+03 


C. 50245+01 


21 


0. 542AE+01 


0.84C0E+02 


0.50245+01 


17 


0. 5A2AE+C1 


C. 1663E+03 


C. 50245+01 


20 


0. 542AE+01 


0.1232E+03 


C.5024E+01 


19 


C.542^E+C1 


0.1232E+03 


0.5024E+01 


2 1 


C. 5424E+C1 


0.6456E+02 


C. 50246+01 


19 


0.542AE+01 


0.64566+02 


0 .50245+01 


2 C 


C. 5424E + 01 


0. 84C0E+02 


0. 5024E+01 


29 


0. E<;745 + C2 


0. lOeiE+03 


C. 19CCE+02 


28 


0.5979E+02 


0.1031 £+03 


C.19C0E+02 


26 


0.597CE+C2 


0. 1112E+03 


C. 19005+02 


27 


0.597SE+C2 


0. 11125+03 


C.19CC5+02 


29 


C.10«8E+C2 


0.21125+03 


C. 99525+01 


2 6 


0. 1C«8E+C2 


C. 21125+03 


C. 99525+01 


32 


0. 5424E+01 


0.12095+03 


C. 50245+01 


3C 


C.5424E+C1 


C. 1209E+03 


0.50245+01 


31 


0 . 109eE+C2 


0.16 26 E+0 3 


C. 99525+01 


31 


C.1098E+02 


0.1626E+03 


C. 99525+01 


24 


0.1098E+02 


0.68^65+02 


€.99525+01 


24 


C. 1G<8E+C2 


0.6846 E+02 


0.99525+01 


26 


0. 1098E+C2 


C. 6846E+02 


C. 99526+01 


29 


0. 1098E+02 


0.6846 E+02 


C.9952E+01 


25 


C. 1C98E+02 


0. 1384E+03 


0.99525+01 


2 5 


0.109eE+02 


0. 1384E+03 


C. 99525+01 


29 


0. 1648E+02 


0.1558E+03 


0.10825+02 


28 


0. 164£E+C2 


G. 15585+03 


€.10825+02 


32 


0. 1648E+02 


0.1240E+03 


€.10825+02 


3C 


C. 1648E+02 


0.12405+03 


0.1082E-*'02 


32 


0. 164€E + C2 


0.76845+02 


0. 10825+02 


31 


0. 1648E+02 


0.768^6+02 


0.10325+02 


29 


0. 164£E + C2 


0.8400E+02 


0.10825+02 


6C 


0.E142E-C6 


0. lOCOE+00 


C.O 


6C 


0.3142E-06 


O.lOOOE+00 


0.0 


32 


C. 16« 6E+C2 


C.96C0E+02 


C. 10825+02 


39 


0. 5979E+C2 


0.96 115+02 


C. 19C0E + 02 


37 


0. 59T9e+C2 


0.96115+02 


0 .19005+02 


34 


C. 5979E+C2 


0.9379E+02 


€. 19005+02 


36 


0.5979E+C2 


0.9879 6+02 


C. 19005+02 


39 


0. Icfc 7E+C2 


0.1792E+03 


0. 10565+02 


37 


C. 1667E+C2 


0. 1792E+03 


C. 10565+02 


^0 


0. 1667E+02 


0.1000 E+03 


0.10565+02 


38 


G. 1667E + 02 


C. lOCOE+03 


€.10565+02 


52 


0. 1667E+02 


0.20665+03 


C. 10565 + 02 


53 


0. 1667E+C2 


C.2066E+03 


C. 10565+02 


9 ^ 


0. 1667E+C2 


0.63S2E+02 


C. 10565+02 


33 


0. 1667E+02 


0.6332E+02 


0. 10565 + 02 


37 


0. 1C67E+C2 


0.63825+02 


0. 10565+02 


- c 


0.1667E+C2 


0. 6382E+02 


C. 1056E+02 


39 


0. 1667E+02 


0.1289 E+03 


0.10565+02 


37 


C. 1667E+C2 


C. 1289E+03 


C. 10565+02 


^0 


0. 1667E+02 


0.32545+C2 


C. 10565+02 


33 


0. 1667E+C2 


0.32545+02 


0.1056S+02 


47 


C. 5979E+C2 


0.72C9E+02 


C. 19005+02 


45 


0.5979E+02 


0.7209E+02 


0.19005+02 


^2 


0. 5979E+C2 


0.7411E+02 


0. 19005+02 


^3 


C. 'C79E+C2 


0. 74115+02 


C. 19005+02 


47 


0. 1697E+02 


0.1446 E+03 


0.10355+02 


4 5 


C. 1697E+C2 


0. 1446E+03 


C. 10355+02 


48 


0. 1667E+02 


0. 7428E+02 


€.10565+02 


46 


0. 1667E+02 


0.74285+02 


0.10565+02 


26 


0. 16675+C2 


C. 71 76E+02 


C. 10565 + 02 


27 


0. 1667E+02 


0.71765+02 


C. 10565+02 


- c 


C. 1667E + C2 


0. 12175+03 


0. 10565+02 


3 c 


0. 1667E+02 


0. 1217E+03 


C. 1056E+02 


35 


0. 1667E+02 


0.61605+02 


0. 10565 + 02 


36 


C. 16675 + 02 


0.6156E+02 


C. 10565+02 


42 


0. 1667E + 02 


0.91385+02 


0. 10565+02 


43 


0. 1667E+02 


C. 91385+02 


0.10 565+02 



0.3016E+00 

0.3C16E+00 

0.3016E+00 

0.3437E+00 

C.3437E-I-00 

0.3437E+00 

0.5437E+00 

0.3437E+00 

0.3437E+00 

0.3437E+00 

0.3437E+00 

0.3437E^K)0 

C.10C25-^01 

0.1002E+01 

0.1002E>01 

0.1CC2F+01 

0.3511E+00 

0.3511E+00 

0.3437E+00 

0.3437E+00 

0.3511E+00 

0.3511E+00 

0.3511E+00 

0.3511E+00 

0.3511E+00 

0.3511E+00 

0.3511E+00 

C.3511E+00 

0.4848E+00 

0.4848E+00 

0.4846E+00 

0.43485+00 

0.4848E+00 

0.4848E+00 

0.4848E+00 

0.0 

0.0 

0.4848E+00 
C.10C2E+01 
G. 10025+01 
0. 10025+01 
C.10C25+01 
0.50265+00 
0. 5026E+00 
0.5026E+00 
0.50265+00 
0. 5C26E+00 
0.5026E+00 
0.5026E+00 
0. 5026E+00 
0.5026E+00 
0.5026E+00 
0.5C26E+00 
0.5026E+00 
0.5026E+00 
0 .50265+00 
0. 1002E+01 
0.10C25+01 
0.10025+01 
C.1CC2E+01 
0.521SE+00 
0.52195+00 
0.5C266+00 
0.50265+00 
0.5026E+00 
0. 50266+00 
0.50266+00 
0. 5C26E+00 
0 • 5026E+00 
0.50266+00 
0. 50266+00 
0.5 0265+00 
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VII . CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

An existing finite element code was expanded to encom- 
pass the more general case of frame members; i.e., six 
degrees of freedom per joint. Combined truss and frame 
structures were designed for minimum weight with multiple 
load conditions considered. 

The displacement method was used for static analysis 
and the subspace iteration method was used for eigenvalues. 

Several examples were considered. In every case the 
code worked as an analysis tool, and significant weight 
reductions were obtained with the coupled optimizer COMMIN . 

The SADX design code has been shown to be a useful tool 
for ship mast optimum design. 

3 . RECOMMENDATIONS 

The following recommendations may be of value for future 
work . 

1. The routines necessary to calculate gradients analyt- 
ically should be added to the code. 

2 . The code should be extended to include plate and 
shear ele.ments and a library of member cross-secrions . 

3. An out of core equation solver should be added. 
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4. The ability to specify multipliers for applying iner- 
tial loads along the three coordinate axes should be added 
in a fashion similar to that used for applying structure's 
own weight as loads. Such an addition would simplify dynamic 
load analysis and design. 



5 . The method of gradient calculation should be depen- 
dent on specific gradients required [Ref. 7] and [Ref. 3]. 

5 . Gradients of frequency constraints would benefit from 
a more efficient algorithm [Ref. S]. 



7. The need for a large scale public str 
tion code still exists . 

3 . The code should be modified to allow 
size design followed by simultaneous optim.um 
optimum geometry design. 



uctural 



optimum 
ms mm s m s 



omtimiza- 



rp 0 0 "vn 



1Z6 



and 
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APPENDIX A 



DATA FILES 



A. INTRODUCTION 

This appendix contains the data files used to create 
the test cases in Chapter VI. Additionally the data file 
for the USER'S guide in complete form is presented. 
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TABLE X 



DATA FILE TRUSS-BRACED CANTILEVER BE/^aM 



SBLCCK A TITLE 1 CARO FORMAT 80H 

1RLSS-BRAC5D CANTILEVER BEAM: EXAMPLE i ^ ^ 



I0LCCK 8 



N68 



NEU8C 



NLC 



CCNTFCL 

NEF 

2 



PARAMETERS 
NJ 



L8LCK 



NCSFLC 



NFREQ 



NSTRES 



NCJ 



NFHASS 



NSTW 



CAROS 

NMT 



NEIG 



FORMAT 3110 
lOVCLC 



NFMW 



NEIGl 



NOJ 



NPRI 



SBLGCK C DYNAMIC ANALYSIS INFORMATION 1 CARO FORMAT I10,2F10.0 
S LMASS . GRAY EPScIG 

5 I I 

0 386 .4 0 . 

♦ 4c4c«5 3»)»c «6 ^ ^ ^ 7 ^ 

$81CCK C JOINT COORDINATES 



NJ 

JN 



S8LCCK E 
< NOJ 
S JN 
$ 



CAPOS 
^ X-CCCRO 

1 
? 

3 
A 
5 

OESIGN 
CAROS 

I 

4 

5 



I 

0 . 

100 . 

200 . 

0 . 
0 . 



Y-COGRO 

I 

0 . 

0 . 

0 . 

150. 

0 . 



VARIABLE LINKING 
lY 

i 1 



Z-CCORO 

0. 

0 . 

0 . 

0 . 

50. 

DATA 



FORMAT I10,3F10 



IZ 



PCX 



FORMAT 
PC Y 

I 

1.0 
1.0 



4I10t3F10.0 

I I 
1.0 1.0 
1.0 1.0 



*2*’»*<=^=****=^2 «5 3|t4: :k:«c:«c :gi 7 



58L0CK 
s E 
$ 



hater ial prcperti es 

RFC SIGMIN 

1.06+^ O.l -2500ol 

2.9E+7 0.3 -36000. 



NMT 
SIGMAX 

25COol 
36 000 



CAPOS 
K EULER 



FORMAT 6F10.0 
POISSM 



I 

2 . 

2 . 



.27 

.27 



J:»t»<c:it:|t3>t2>l2j:<c3f: * 5 ;(t ^ ^ 5 ;jc:jc :#t *:*:*:*:<:* 7 ^>: 1 : :* ^ *:(t 



S8LCCK G 
S LNO 

5 



EAR ELEMENT 
NCCE2 

1 1 

2 2 



INFORMATION NEB CARDS FORMAT 

NC0E3 MATCOD NSOGl AREA 

i I I 

4 11 

5 12 



51 lo.FiOt no 

I I 

0 1 

3.0 1 



4 :t 2 4 3 4i4:^4c4esc 4c4c4c44(« 4^^ t: 5 4c «4c ^4(*4c 4: j>4c 74 ^ 4 :^#^'^ 



S8LCCK H FRAME ELEMEM INFORMAL 



LNO 



N0CE2 



NCDE3 



I 



CHAPCIH2 

1 

1«0 

2 

1 .0 



ION 2*NEF CAROS FQR'M 

HATCOD NSDGl NSDG2 

1 I 

2 3 



7110/ 2F 10 

I 1 



JCFAROIMl 
S I 

5.0 

4 

5.0 

5* 4:4:3»^:> 4 :414c4c4i ^:44c:»4c 4t 2 44 t 4: 4t * * 4 c3 4c4c4c3>^4c 4c:># 4*:* *4c 4c:» :>4:4c 5 4:4t4: 4« 4.'4c 4: 4c6 -T4cC 4c4c 4t C4c 7 4t 4c ^ ^4i3^ 4 4c 4c 



58 LOCK 
S JN 
5 



I JOINT CC2^STPAINT DATA 
IX lY 

i 1 



IZ 



NCJ 

I 

I 
0 
0 
1 
1 



CAROS 
IXX 



FORMAT 
I YY 



7110 



0 

0 

1 

1 



IZZ 



0 

0 

1 

1 



4 44c4c*4:ti3*:** 4 4*404 41 4*4c4c ** **4c A** 44: 344: <c 5 4t* 4c :4 ** * * * 6 <^=» 7 * ** 

SBLGCK J JCIMT LOACI NG CATA 1 CAPO FORMAT 110 

i (CMIT IF NLC=C IK 8LCCK B» 

$ NLJ 
S I 

1 
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TABLE XI 



DATA FILE TRUSS-BRACED CANTILEVER BEAM continued 



I 

0 • 

- 1000.0 



NLJ CAROS 
TX 



I JN I FX j '"'^1 

3 0. lOCO.O 

3 0. 0.0 

SBLGCK K FIXED MASS CATA NFMASS CA 

5 (CMT IF KFHASS = 0 IN BLOCK 8) 

t JN MASS 

$ I I 

3 250. C 

I* ^ ^ 4 * J 3|t ^ ^ 5 * * * ** * 4c 4** 

SBLCCK L DESIGN VARI>5PLF I.NPORMATION CARO 



FORMATIXO,6FlO 

TY TZ 

J I I 

0 * 0 . 0 . 



0 « 0 « 0 « 

4 A # 5 4( 4e 9 ^ ^ ^ 4c 4 « 7 » « * « 4( ^ ^ « 

ROS «=ORMAT IlOtFlO.O 



(AREA AND CIMENSICN VARIABLES 



CMT 



:><: 4cC*iCe4c^:*64'*^4c 

S AS REQn FORMAT 3F10.0 



$ 

s 


XA( n • 


. . 

1 


XA(NOVARlj 


1 


1 




2C.C 


20.0 


20.0 


2.0 


2.0 


$ 

s 


XAL( I) 


1 


1 


1 


1 




.50 


. 50 


4.0 


.10 


• 10 


$ 

s 


XAU( I| 
25.0 


1 

35. C 


1 

25.0 


2.5 


1 

4.0 



3 :»4t4c4c 44 c 5 #4t ^ 4c ^ 

SBLCCK M DESIGN VARIABLE INFORMATION CAROS AS REO*D FORMAT 8F10.0 



(CCOROINATE variables 



XC(I 



f • 





150.0 


50.0 


$ 


XCL(I) 




5 


1 

0.0 


1 

0.0 


S 


XCLMI) 




s 


I 

200.0 


1 

ICO .0 



I XC(N0VAR2j 



I 



QMI 

I 

I 

I 



t# «**♦♦** 1 ««::»*<*»* *2 **=»=*-■**» =*3 
561CCK N JOINT C!S?LAC=MENT CONSTRAINTS 
S (CMIT IF NOSPLC=C IN BLOCK 3) 

t JN . CIR , LC . SL 

5 I I I I 

3 2 1-3.0 

3 3 2 -3. 5 

5»*»***4«1 ■^7^*9»**»2*^^***^'*^3 
3BLCCK C FOEOUSNCY CONSTRAINTS NO. OF C 
! (CM IT IF NFRSC=0 IN BLOCK 3) 

I PC^ . .PCN 

1.0 

<:> 4c# 4: 4 =>* 4 l# # 2 4t :>3 4*4c# *4t A*’* 

SELCCK F EMO CARD 
END 



T IF NOVAR2=OI 

I 

I 



44 C 4 C # 4 t 4 4 c 54 '*=> 4 t 4 4 s 4 c 4 46 =>*^ 4 t # 4 t # 4 c 4 ( 74 ^: 44 ! 4 c 4 c 4 c ♦♦♦ 

NOSPLC CAROS FORMAT 3110, 2F10 



BU 



I 



3.0 

3.5 

# 4 c 4 c # 4 c # 4 c 5 4 :^ 4 c 4 c # 4 c 4 c 4 c 4 ^ 4 c 4 c 4 c 4 c # 4 : # 4 c * 7 ** * 4 c 4 c 4 c 4 f 4 ! 4 c 

AP.OS AS RE0*0 FORMAT 8F10.0 



I I I 

»««»#* 5 *»:#*»**** 6 ****< t :* ac #<[ 7 ********* 
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TABLE XII 



DATA FILE TWO-TIER 3-D PORTAL FRAME 



52 LCCK A title 1 CARO FORMAT 80 H 

2 -TIEF FRAMJ= STRUCTURE WITH TRUSS X-BRACcS: EXAMPLE 2 

5 ♦★*44:^ *4i«3* *5*4*5 4i4: *****6*********7****** *** 



5 ELGCK e CONTROL PARAMETERS 



NEB 



NEUBC 



NLC 



I 

20 



NEF 



L 31 .CK 



16 



NJ 



NFREQ 






NOSFLC 



li 



NSTRES 



NCJ 



NFMASS 



NSTW 



3 CAROS 
I NMT 
8 

NEIG 



FORMAT 8110 

lOVCLC NOJ 



NEIGl NPRI 



I 



NFMW 



|********X*** ****** 2 *********4** *******5 *********£)********* 7 ****** *** 
SELCCK C CYNAMIC ANALYSIS INFORMATION 1 CARO FORMAT I 1 C, 2 F 10.0 
S LMASS GR/V EPSEIG 

$ I I I 

0 3 66. A 0. 

t* **** ***!** ****** *2 *’>**♦<'*^*3 *:<t*4t*****4*********5*********6*********7********* 

S 3 L 0 CK C JOINT CCOROINATES 



NJ CAROS 
JN X-CCQRO 



Y-COORO Z-COORO 
1 



FORMAT I 10 » 3 F 10 



1 


100. 


0. 


100 


2 


100. 


0. 


•100 


3 


-100# 


0. 


•100 


4 


•100. 


0. 


100 


5 


100. 


50. 


100 


6 


100. 


50. 


• 100 


7 


•100. 


50. 


•100 


8 


• 10 c. 


5 0. 


100 


9 


100. 


100. 


100 


10 


10 0. 


100. 


•100 


11 


•100. 


100. 


•100 


12 


•100. 


100. 


100 



I***** *** X«* ^*** 5 ?s<: 2***’>**^*^3 *** 4 ** **** *** 5 ** **** ***^ ********* 7 ** ♦*♦* *♦* 

SELCCK E CESIGN VARIABLE LINKING DATA 

S NCJ CARDS FORMAT AI 10 , 3 F 10.0 

J JN I> . lY . IZ , PCX , PCY , PCZ , 

J ] I I I • , i 

1 1 0 2 1.0 0.0 1.0 

2 1 0 2 1.0 0.0 - 1.0 

3 10 2 - 1.0 C.O - 1.0 

^102 - 1.0 0.0 1.0 

5 3 5 1.0 1.0 1.0 

6 3 4 5 1.0 1.0 - 1.0 

7345 - 1.0 1.0 - 1.0 

8 3 4 5 - 1.0 1.0 1.0 

** *»:****»2 >**3»*****3 ****:t»:***4*^»*** **»5 at* # 3 *^,»*** *»«»* 7 *»»♦*♦ *»» 

IBLOCK F MAT 5 RIAL PRCP 6 RTIHS NM 7 CAPOS FORMAT 6 F 10.0 

S E RFC SIGMIN . SIGMAX . K EULER . POI SSN 



2.5E+!r 0.3 -42COol 36COoi 4 I .2^ 

5 ********X*******<s^ 2’*^^*^*^^^3 *********4** **** ***5**** ***** 6 *** ****** 7********* 
ISLGCK G BAR ELEMENT INFORMATION NEB CAROS FORMAT 5110. FIO. 110 



LNC 



I 



N 0 CE 2 



2. 5. 2. 1.1. 5. 0. 1 

3.2.7 .1 ,1 ,5.0,1 

4. 6. 3. 1.1. 5. 0. 1 

5. 3.8.1 ,1,5.0, 1 
e, 7, 4,1 ,1,5.0, 1 
7,4,5, 1, 1,5. 0, 1 

6. 1.8. 1.1. 5.0. 1 

5. 5.7. 1. 1.5. 0, 1 

10 .6.8. 1.1. 5. 0. 1 
11 , 5 , 1 c ,1,2,5. C,1 
12 ,6,5, 1,2,5. C,1 
13,6, 11, 1,2, 5. C,1 
14, 7, 1C, 1,2, 5. C,1 



1 



NCDE3 



MATCOC 



i 



NSDGl 



AREA LSECT 

5.0 
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TABLE XIII 



DATA FILE TWO-TIER 3-D PORTAL FRAME continued 



15 , 7 , 12 , 1 , 2 , 5 . C,l 

16 . 8 . 11 . 1 . 2 . 5 . C ,1 
17 , 6,9 , 1 , 2 , 5.0 a 

18 . 5 . 12 . 1 . 2 . 5 . C ,1 
19 , 9 , 11 , 1 , 2 , 5 . 0,1 
20 fl 0 vl 2 tl ?2 ?5 •Ofl 

SBLCCK h PPA.VE ELEMENT INFORMATION 
S LNC NCC 62 NCOE 3 MATCOC 



I 



i 1 I 

21 ,lt 5 ,l. 3 ,A,l 

5 .C 0.50 

22f2»6»lf3f4, 1 

5.0 0.50 

23 » 3 , 7 ,lf 3 » 4 tl 

5.0 0.50 

24 t 4 , 8 tlf 3 t 4 ,l 

5 .C 0.50 

25 f 5 , 9 ,l, 5 , 6 ,l 

5.0 0.50 

26,6tlCfl»5fEf 1 

5.0 0.50 

27 , 7,11 tl» 5 » 6 fl 

5.0 C .50 

28 . 8 . 12 . 1 . 5 . 6.1 

5.0 • 0.50 

29 . 5 . 6.1 , 7 , 8, 1 

5.0 0.50 

30 , 6 , 7 ,1 , 7 , 3,1 

5 .C C .50 

31 . 7 . 3 . 1 . 7 . 8.1 

5.0 0.50 

32 . 6 . 5 . 1 . 7 . 8.1 

5.0 • 0.50 

33 . 9 . 10 . 1 . 9 . 10.1 

5.0 C .50 

34 . 10 . 11 . 1 . 9 . 10.1 
5.0 0.50 

35 . 11 . 12 . 1 . 9 . 10.1 
5.0 0.50 

36 . 12 . 9 . 1 . 9 . 10.1 
5 .C 0.50 



5 BLGCK I JOINT CCNSTPAINT DATA 
J JN . I> . lY 

S I I 

1 , I, 1 , 1 , 1 , 1 , 1 

2 , 1 , 1 , 1 , 1 , 1, 1 

3 , 1 , 1 , 1 , 1 , 1,1 

Afl, 1,1, 1,1,1 

5 , 

6 , 

7 , 

6 , 

9 , 

10 , 

11 , 

12 , 



SBLCCK J JOINT LCACING CATA 
% (CMT IF NLC =0 IN BLOCK B) 

I NLJ 



2 ^NEF CAROS «=ORMAT 7 I 10 / 2 F 10 
NSDGl I NSDG 2 ^ LSECT ^ 



IZ 



NCJ CARDS 
I XX 

I 



FORMAT 7110 

I I I 



s 
2 , 
i 

S JN 
S 

9 , 1000 . 

10,1000. 

^ , 

9 . 500 . 

10 . 500 . 



I 



FX 



FY 



FZ 



1 CARD 



NLJ CAPOS 
TX 
I 



FORMAT no 



FORMAT I 10 , 6 F 10 
TY TZ 

I 
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TABLE XIV 



DATA FILE H^JO-T^ER 3 -D PORTA.L rRAI-lE continued 



5. -5CO. 

6, -SCO. 

1C tiooo. 

<:.C.CfC.CflOOO. 

12 .-lOGC. 
11^C.0»0.C,-10C0. 

I I 



JN 



1 
1 

5C.C 
lOtSO.C 
11 ,50.0 
12 ,50. C 

58LCCK L 



CESIG^ VARU3LE INPCR'^ATICN 



(AREA ANC CIMENSICN VARIAELES 



XA( I| 

3.5 
5.0 
XAU I| 

G.l 
.50 
XAU( I j 

3C.0 
3C.C 



3 

0 . 



J 



;«A(NOVAR Ij 
5.0 



CAROS AS RS'3'D FORMAT 8F10.0 
CHIT IF NDVAR1 = 0) 



oA 



0.1 


.50 


0. 1 


0.1 






1 


1 


i 


30.0 


30 .0 


.50 



setOCK y OESIG^ VARIABLE INPGRMATICN 



(CCCR0INAT5 VARIABLES 



XC( I j 



I 



XC(N0VAR2} 



CARDS AS 
3MIT IF 



s.i 

.si 

I 

30.0 



RcQ'O 

N0VAR2=0) 



oA 



0 . 1 

I 

.50 
FORMAT 



5.i 



I 

.50 



I 

30. C 



8F10 .0 



C.5 

O.l 

c.s 



ICC.O 
XCLd j 

20. 0 

XCU(i| 

16C.0 

SBLCCK N JOINT DISPLACEMeNT 



ICO.C 


ICO.O 


ICO.O 


50.0 


20 A 


I 

20.0 


1 

20.0 


1 

20.0 


160 .<i 


160.0 


160. i 


1 

30.0 






t ^ * 4:4 zfc a ^ 


aaac5 aa 



IF NCSPLC=0 IN 3L0CK 3 ) 



S (CMT _ 

1 JN . OIR LC 

i 1 

9 

10 
11 
12 
9 

10 
11 
12 
9 
10 
11 
12 

SELCCK C FPcCUENCY CONSTRAINTS 
$ (CyiT IF NFREC-0 IN BLCCK 

5 FCl . . FCN 

Sill 
l.C 5.0 10.0 

5ELCCK f END CARC 
ENC 



:nstra:nts ndsplc cards format 5rio,2Fio 



1 
1 
1 
1 
1 

2 
2 
2 
2 
3 
3 
3 
3 



3L 



I 

-5.0 
-5.0 
-5. 0 
-5. 0 
-5.0 
-5.0 
-5.0 
-5.0 
-5.0 
-5.0 
-5.0 
-5.0 



BU 



5.0 

5.0 

5.0 

5.0 

5.0 

5.0 

5.0 

5.0 

5.0 

5.0 

5.0 

5.0 



3i 



NO. OF CAROS AS REQ’D 



«=GRMAT 3F10.0 



P 



TABLE XV 



DATA FILE DD-963 FOREMAST 



JBLOCK A TITLE 

CC-963 CLASS CESTROYEP FCPHARO MAST REDESIGN: EXAMPLE 1 

$ THIS IS A SI^PLIFlCiTIC^ OF THE ACTUAL STRUCTURE Af40 LOADS BASED ON 
i NAVSMPS 0P4V.IN6 NUMBER 128-A53S510 AND N'AVSHIPS SKETCH NUMBER 
1 8006A-128-SH4537<; 79 

4 THIS ANALYSIS RECLIPES A THE X-AXIS FORWARD, 

J TO ST60; A CIFFERENT ORIENTATION THAN IN THE 



NEB 



5ELCCK e CCNTSCL PARAMETERS 
NEF 

172 
UBCCK I 

1 

NOSFLC 



Y-AX IS UP, AND Z-AX IS 
ABOVE DRAWINGS 



i NEUBC 



NLC 



li 



NJ 


NCJ 


NMT 


73 


1 

20 




NFREQ 

1 


NFMASS 

i 


NEIG 


0 


33 




NSTRES 

1 


NSTW 

i 


^F^(W 



1 



I 

1 

I 

0 

I 

1 



IDVCLC 



NEIGl 



NDJ 



NPRI 



SELCCK C CYNANIC ANALYSIS INFORMATION 
C 386.4 0. 

SBLOCK C JOINT CCOROINA7ES 

-1C6.0 
1C6.0 
42.0 
-42.0 



-240.0 
-240.0 
O.C 
0 .0 



0.0 

0.0 

0.0 

0.0 



$ NOCE 



5 


-223. C 


54.0 


0.0 


6 


-120. C 


54.0 


68.4 


7 


-2.64 


54.0 


0.0 


8 


-120.C 


54.0 


-68.4 


9 


-216. C 


1C8.0 


-94. 7 


10 


-216 .0 


103 .0 


94.7 


11 


-5. 28 


1C8.0 


42. 0 


12 


-5. 28 


1C8.0 


-A2.0 


13 


-7.92 


162.0 


0.0 


14 


-172. C 


444.0 


0. 0 


5 IS 


AN UNCCNNECTEC FIXED 


DUMMY N 


15 


0.1 


0.0 


0.0 


16 


-192.6 


216.0 


-83.16 


17 


-192.6 


21fi> .0 


83.16 


18 


-1C. 56 


216.0 


42.0 


19 


54.0 


216.0 


42.0 


20 


54.0 


216.0 


-42.0 


21 


-1C. 56 


216.0 


-42. 0 


22 


-172 .C 


736.0 


- 18.0 


23 


-172 .0 


736.0 


18.0 


24 


-13.2 


270.0 


0. 0 


25 


-168 .8 


324.0 


0.0 


26 


-163. E 


524.0 


-71.76 


27 


-168 .8 


324.0 


71.76 


28 


-15. 84 


324.0 


42.0 


29 


-15. 84 


324.0 


-42.0 


30 


1C8 .0 


324.0 


48.0 


31 


168 .0 


324.0 


0.0 


22 


1C8.C 


324.0 


-43. 0 


33 


-18. 12 


372.0 


0.0 


34 


-147 .8 


420.0 


-6 1.6 8 


^ c 


-147.8 


420.0 


0. 0 


36 


-147.8 


420 .0 


61.63 


21 


-20.4 


420.0 


42.0 


“^8 


12 .C 


420.0 


39.0 


39 


-20.4 


420 .0 


-42.0 


4C 


12.0 


420.0 


-39.0 


41 


-22 .2 


456.0 


0.0 


42 


-132 .0 


492.0 


-54.0 


42 


-132. C 


492.0 


54. 0 


44 


-172 .0 


492.0 


- 54.0 


45 


-24. C 


492.0 


42.0 


46 


20.0 


492.0 


36. 0 


47 


-24.0 


492.0 


-42.0 


48 


30. C 


492.0 


-36.0 



$ END NCC5S FCP THE LOWER QUAORAPOD CF THE 
S START THE UPPER FCFEKAST NGOES 



FORWARD MAST 
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TABLE XVI 





DATA 


FILE DD- 


96 3 FOl 


49 


-172. C 


492.0 


-120.0 


50 


-172. C 


492.0 


120.0 


51 


-228.0 


408.0 


0.0 


52 


-288.0 


492.0 


-56.0 


53 


-288 .C 


492.0 


56.0 


54 


-240.0 


492.0 


0.0 


c c 


-172. C 


492.0 


-2 76.0 


56 


-172 .0 


492.0 


276.0 


57 


-213 .0 


600.0 


-24.0 


58 


-218.0 


6CO.O 


24.0 


59 


-172 .0 


492 .0 


54.0 



S NCOE IS 4N L'NCCNNECTEC FIXED CU^MY MODE 



60 


0.0 


0.0 


0.0 


61 


-172.0 


600.0 


0.0 


62 


-172.0 


492.0 


0.0 


63 


- 1 72 . C 


664.0 


0.0 


64 


-172 .0 


736.0 


0.0 


65 


-172. C 


896.0 


0. 0 


66 


-172 .0 


976.0 


0.0 


67 


-172.0 


1056.0 


0.0 


68 


- 1 72 . C 


736.0 


-2 64. 0 


69 


-172 .0 


736.0 


264.0 


7C 


-172 .C 


736.0 


-156.0 


71 


-172 .0 


736.0 


156. 0 


72 


-172 .0 


736.0 


-228.0 


72 


- 1 72 . C 


736.0 


228. 0 



4SIQCK F yATEF 
2.S5+7 

5# ^ 4 1 4> 

IBLGCK G BAR E 
$ LNC UC 

S I 

5BLCCK h FRAME 
1 START FRAME 
1 

11.75 

2 

11.75 

a 

11.75 

4 

11.75 

c 



4 ^*^^^**5^^^^*4i^*e****^**^^ ?*♦♦ 

lAt INFCPyATICN 

0.1 -A2000. 42C00. 4 

LEM6NT INFCRMATIGM NEB CARDS 

C£2 NC0E3 MATCOC NSCGl 

I I t 

4 ^ 4t 2 4 4 4 k 4^:$: it 3 « aMt 4r 4t 4i 44c « ^ « 4c3|c 4c * 

ELEMENT INFORMATION 

ELEMENTS FCR 1ST CR LOWEST TIER OF QUAORAPOC 
4 

.75 



.27 

FORMAT 5110, FlOt 110 
AREA LSECT 

* * I 

5 4:4k aK4e 4 k « ^ ^ t * 



3 
. 75 

1 

.75 

2 

.75 

4 



12 

11 

9 

10 

3 



7.625 


.375 






6 


3 


6 


1 


7.625 


.375 






7 


1 


3 


1 


7.625 


.3 75 






8 


2 


6 


1 


7.625 


.375 






9 


8 


12 


1 


7.625 


.375 


11 




10 


6 


1 


7.625 


.375 






11 


8 


9 


1 


7.625 


.3 75 






12 


6 


10 


1 


7.625 


.375 






13 


3 


7 


1 


7.625 


.375 






14 


4 


7 


1 


7.625 


.375 


11 




15 

7.625 


7 


1 


.375 


12 




16 


7 


1 


7.6:5 


.3 75 






17 


1 


5 


1 


7.625 


.3 75 






18 


2 


5 


1 


7.625 


.375 







1 

1 

1 

1 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 



2 

2 

2 

2 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 



33 



DATA 



T-IZ DD-96 3 JORZMA.ST continued 



19 

7.625 


5 

.375 


9 


1 


3 


4 


2C 

7.625 


c 

.3 75 


10 


1 


3 


4 


21 

7.625 


12 

.375 


9 


1 


3 


4 


22 

7.625 


11 

.375 


10 


1 


3 


4 


22 

7.625 


11 

.375 


12 


1 


3 


4 


2A 

7.625 

ENC LET TIER 


9 10 

.375 

CF CU;iCP4?0C START 


1 

2NC TIER 


3 


4 


25 

11.75 


12 
. 75 


21 


1 


1 


2 


26 

11.75 


11 
. 75 


18 


1 


1 


2 


27 

11.75 


9 

.75 


16 


1 


1 


2 


28 

11.75 


1C 
. 75 


17 


1 


1 


2 


29 

7.625 


9 

.375 


21 


1 


5 


6 


30 

7.625 


10 
.3 75 


18 


1 


5 


6 


31 

7.625 


11 
.3 75 


13 


1 


5 


6 


22 

7.625 


12 

.375 


13 


1 


5 


6 


33 

7.625 


13 

.375 


18 


1 


5 


6 


34 

7.625 


13 

.375 


21 


1 


5 


6 


t c 

7.6^ 


1C 

.375 


16 


1 


5 


6 


36 

7.625 


21 

.375 


16 


1 


7 


3 


37 

7.625 


17 

.375 


13 


1 


7 


8 


28 

7.625 


18 

.375 


21 


1 


7 


3 


29 

7.625 


16 

.375 


17 


1 


7 


a 


40 

11 .625 


12 
.3 75 


20 


1 


7 


9 


11.625 


11 

.375 


19 


1 


7 


3 


^2 

7.625 


20 

.375 


21 


1 


7 


3 


43 

7.625 


13 

.375 


19 


1 


7 


3 


44 

7.625 


19 

.375 


20 


1 


7 


8 


ENC ThE 2ND 


TIER C? ThE 


CUAORAPOC START THE 


3R0 TIER 




45 

11.75 


21 

.75 


29 


1 


1 


2 


46 

11.75 


18 
. 75 


28 


1 


1 


2 


47 

11.75 


^ 4 

.75 


26 


1 


1 


2 


43 

11.75 


17 
. 75 


27 


1 


1 


2 


4C 

7.625 


16 
.3 75 


29 


1 


9 


10 


5C 

7.625 


17 

.375 


23 


1 


9 


10 


51 

11 .625 


20 
.3 75 


32 


1 


7 


3 


52 

11.625 


19 

.375 


30 


1 


7 


3 


7.6 25 


20 

.375 


31 


1. 


9 


10 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

l 

1 

1 

l 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

l 

1 

1 

1 

l 

1 

1 

1 



TABLE XVIII 





DATA FILE 


DD- 


963 


FOREMAST 


continued 




54 


19 




31 




1 


9 


10 


1 


7.625 


.3 75 














55 


18 




24 




1 


9 


10 


1 


7.625 


.3 75 














56 


21 




24 




1 


9 


10 


1 




.375 










57 


24 




28 




1 


9 


10 


1 


7.625 


.375 








5a 


24 




29 




1 


9 


10 


1 


7.625 


.3 75 












59 


16 




25 




1 


9 


10 


1 


11.50 


. 5C 














60 


17 




25 




1 


9 


10 


1 


11. 5C 


. 5C 
















MEMBERS 61 


ThRU 70 


MCCEL 


THE 


SPQ- 


9 PLATFORM 


ALONG WITH 


A FIXED 


MASS 


61 


26 




29 




1 


11 


12 


1 


11. 5C 


. 50 














62 


27 




28 




1 


11 


12 


1 


11.50 


.50 
















£3 


29 




32 




1 


11 


12 


1 


11.50 


.50 
















64 


20 




30 




1 


11 


12 


1 


11.50 


. 50 
















65 


31 




32 




1 


11 


12 


1 


11.50 


. 50 
















66 


31 




30 




1 


11 


12 


1 


11.50 


.50 
















67 


28 




29 




1 


11 


12 


1 


11.50 


.50 
















ELEMENTS 63 


/NO 69 


APE DUMMY 


ELEMENTS 








68 


15 




60 




1 


0 


0 


1 


.CCl 


.0001 
















69 


15 




60 




1 


0 


0 


1 


.CCl 


.0001 
















70 


30 




32 




1 


11 


12 


1 


11. 5C 


. 50 
















ENC Tl-E 3RD 


TIER CF 


THE 


OUAORAPOC 


START THE 


4TH TIER 






71 


29 




39 




1 


1 


2 


1 


11.15 


.75 
















72 


28 




37 




1 


1 


2 


1 


11.75 


.75 
















13 


26 




34 




1 


1 


2 


1 


11.75 


.75 
















74 


27 




36 




1 


1 


2 


1 


11.15 


. 75 












14 




75 


26 




39 




1 


13 


1 


7.625 


.375 














1 


76 


27 




37 




1 


13 


14 


7.625 


.3 75 












14 


1 


17 


29 




40 




1 


13 


7.625 


.3 75 
















70 


28 




38 




1 


13 


14 


1 


7.625 


.375 










13 


14 




79 


26 




52 




1 


1 


7.625 


.375 










13 


14 




eo 


27 




53 




1 


1 


7.625 


.3 75 










13 


14 




61 


28 




33 




1 


1 


7.625 


.3 75 










13 


14 


1 


82 


29 




33 




1 


7.625 


.375 










13 


14 




83 






37 




1 


1 


7.625 


.375 












14 




84 


- -a 




39 




1 


13 


1 


7.625 


.3 75 










13 


14 


1 


65 


34 




39 




1 


7.625 


.3 75 










13 


14 


1 


86 


36 




37 




1 


7.625 


• 375 










13 


14 


1 


87 


39 




40 




1 


7.625 


.375 










13 


14 


1 


88 


37 




38 




1 
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124 
7.625 

125 
7.625 

126 
7.625 

127 

l.tll 

123 

7.625 

i ENC ThE LO»^ER 

129 
7.625 

130 
7.625 

121 

7.625 

122 

7.625 

123 

7.625 

134 

7.625 

125 
7.625 

126 
7.625 

127 
23. CO 

13E 
23. CO 

129 
23. CC 

140 
23. CO 

141 
23. CO 

1^2 

9.50 

143 
9.50 

144 
9.50 

145 
7.625 

! UPPER 'tAROARP 

146 
14. CC 

147 
14. CO 

148 
IC.CC 

149 
IC.CC 

150 
10. CO 

151 
10. CO 

152 
7.625 

153 
7.625 

154 
5.75 

155 
5.75 

156 
7.625 

157 
7.625 

158 



.375 

53 54 

.3 75 

42 54 
.375 

43 54 
.375 

42 44 
.375 

43 59 
.375 

CJADRAPQO START THE 

52 58 

.3 75 

52 57 
.375 

53 58 
.375 

57 61 
.375 

58 61 
.3 75 

57 58 

.3 75 

57 63 
.375 

58 63 
.375 

14 25 

I .5 

14 62 

I .5 

62 61 

1.5 

61 63 

1.5 

63 64 

1 *5 

64 65 
. 50 

65 66 
. 5C 

66 67 
. 50 

27 29 

.3 75 

S ANC BRACES 

64 70 

.75 

64 71 

.75 

70 72 
. 50 

71 73 
. 50 

72 68 
•50 

73 69 
•50 

61 72 

•375 

61 73 

.3 75 

63 70 

• 25 

63 71 

. 25 

61 72 

•3 75 

61 73 

•3 75 

63 70 
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15 


16 


1 


1 


15 


16 


1 


1 


15 


16 


1 


1 


15 


16 


1 


1 


15 


16 


1 
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1 
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17 


18 


1 


1 


17 


13 


1 


1 


17 


18 


1 


• 1 


17 


18 


1 


1 


17 


18 


1 


1 


17 


13 


1 


1 


17 


18 


1 


1 


17 


18 


1 


1 


19 


20 


1 


1 


19 


20 


1 


1 


19 


20 


1 


1 


19 


20 


1 


1 


19 


20 


1 


1 


21 


22 


1 


1 


21 


22 


1 


1 


21 


22 


1 


1 


13 


14 


1 


1 


23 


24 


1 


1 


23 


24 


1 


1 


25 


26 


1 


1 


25 


26 


1 


1 


25 


26 


1 


1 


25 


26 


1 


1 


27 


28 


1 


1 


27 


28 


1 


1 


29 


30 


1 


1 


29 


30 


1 


1 


27 


28 


1 


1 


27 


28 


1 


1 


29 


20 


1 
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T.£25 .375 



: ThF 


KEil CF THE CJiCRiPOC 


START 


THE 5TH TIER 






59 


3 9 47 


1 


1 


2 


1 


11.75 


.75 










<iC 

11.75 


37 45 

. 75 


1 


1 


2 


1 


91 

11.75 


34 42 

. 7 3 


1 


1 


2 


1 


92 


26 43 


1 


1 


2 


1 


11.75 


. 75 








c-a 


34 47 


1 


15 


16 




7.^i5 


.3 75 






94 

l.i2S 


36 45 

.3^5 


1 


15 


16 


1 


cc 


40 43 


1 


13 


14 


1 


7.625 


.375 






9c 


35 46 


1 


13 


14 


1 


7.625 


.3 75 






c7 


25 26 


1 


13 


14 


1 


7.625 


.375 






93 


25 27 


1 


13 


14 


1 


7.625 


.275 






99 


26 35 


1 


13 


14 


1 


7.625 


.375 




ICC 


27 35 


1 


13 


14 


1 


7.625 


.3 75 






ICl 


34 25 


1 


13 


14 


1 


1 • tl^ 


.3 75 






1C2 


35 26 


1 


13 


14 


1 


7.625 


.375 






IC3 


35 42 


1 


13 


14 


1 


7.625 


.375 






1C4 


3 5 43 


1 


13 


14 


1 


7.625 


.3 75 








1C5 


37 41 


1 


13 


14 


1 


7.625 


.375 










1C6 


29 41 


1 


13 


14 


1 


1 milt 


.375 








1C7 


41 45 


1 


13 


14 


1 


7.62* 


.375 








1C3 


41 47 


1 


13 


14 


1 


7.625 


.375 
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-60 
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A 


1/) 

-1 

o 
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X 

II 
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42 47 


1 


1 5 


16 


1 


l.t25 


.375 










lie 


43 45 


1 


15 


16 


1 


7 .625 


.375 










111 


42 43 


1 


15 


16 


1 


7.625 


.375 










112 


45 47 




15 


16 




7.625 


.375 
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4c 47 


1 


15 


16 


1 


7.625 


.375 










11^ 


45 46 


1 


15 


16 


1 


7.625 


.2 75 












43 46 


1 


15 


16 


1 


7.625 


.3 7 5 




15 






116 


52 53 


X. 


16 


1 


7.625 


.375 










117 


26 51 


1 


13 


14 


1 


7.625 


.375 










115 


27 51 


1 


13 


14 




*^.625 


.375 










119 


52 51 


1 


13 


14 


1 


7.625 


.375 










120 


53 51 


1 


13 


14 


1 


7.625 


.375 






16 




121 


4^ 52 


1 


15 


1 


7.625 


.375 




15 






122 


53 59 


1 


16 


1 


7.623 


.3 75 




15 






122 


52 54 


1 


16 
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5.75 


.25 


159 


63 


5.75 


.25 


5 FLATFQFM 




160 


63 


5.75 


.25 


161 


63 


5.75 


.25 


162 


22 


5.75 


. 25 


163 


64 


5.75 


.25 


164 


64 


5.75 


.25 



LOViER ^[AROAR^ 

165 
7.625 

166 
7.625 

167 
11. fO 

168 
11.50 

165 

11.50 

17C 
11. 5C 

171 
5.75 

172 
5.75 



AND ePACES 

.375 
59 
.375 
49 
. 50 
59 
.50 

49 
. 50 

50 
. 5C 

14 

.25 

14 

.25 



71 



22 

23 

23 

23 

22 

62 

62 

44 

50 

55 

56 

49 

50 
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JN 



1 
1 

2 
■a 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 
19 
2C 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

•3 C 

36 

37 






lY 



1 

1 

1 

0 

C 

0 

0 

c 

0 

0 

c 

0 

0 

1 
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0 
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0 

c 

0 

0 

0 

0 

0 

0 

0 

c 

c 

0 

0 

0 

0 

0 

0 

0 



1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 



NCJ CAROS 
I XX 
I 
1 
1 
1 
1 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
1 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 



29 

29 

29 

29 

29 

29 

15 

15 

31 

31 

31 

31 

33 

33 
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I 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



30 

30 

30 

30 

30 

30 

16 

16 

32 

32 

32 

32 

34 

34 

7110 



1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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i 

$ 

« 

s 



5 

$ 

S 



JN 




FX 


FY 


l-ALYARC AND 


hIRE A^TE^NA LCAOS 






-5C. 


-100.0 




50 


-50. 


-100.0 




55 


-50. 


-100.0 




56 


-50. 


-ICO.O 




a 


0. 


-50.0 


SHIP • 


S ^^OTIC^ LOADS 




JN 


1 


FX 

1 


FY 

1 




1 

1 


0.0 


0.0 




2 


o.c 


0.0 




3 


0.0 


0.0 




% 


o.c 


0.0 




5 


0.0 


0.0 




6 


0.0 


0.0 




7 


0.0 


0.0 




8 


0.0 


0.0 




9 


o.c 


0.0 




10 


0.0 


0.0 




11 


0 .0 


0.0 




12 


o.c 


0.0 




13 


0.0 


0.0 




14 


0 .0 


0.0 




16 


0.0 


0.0 




17 


0.0 


0.0 




18 


0.0 


0.0 




19 


0.0 


0.0 



NL J 


CAROS 


FaRMATI 10 . 6^10 






FZ 


TX 




TY 


TZ 


1 

0 . 


1 

0.0 




1 

0 . 


1 

0 . 




0.0 




0 . 


0 . 




0 . 


0.0 




0 . 


0 . 




0 . 


0.0 




0 . 


0 . 




0 . 


0.0 




0 . 


0 . 




c. 


1 


TX 


j 


TY 

1 


TZ 


1 


160.5 




0.0 


0.0 




0.0 


163 . 3 




c.o 


0.0 




0.0 


< 36.4 




0.0 


0.0 




0.0 


137.5 




0.0 


0.0 




0.0 


1 C 8 . 2 




0.0 


0.0 




0.0 


184.3 




0.0 


0.0 




0.0 


66 . 7 




0.0 


0.0 




0.0 


1 64.3 




0.0 


0.0 




c.o 


414.9 




0.0 


0.0 




0.0 


467.0 




0.0 


0.0 




0.0 


381.6 




0.0 


0.0 




0.0 


361.6 




0.0 


0.0 




0.0 


70.2 




c.o 


0.0 




0.0 


325.4 




0.0 


0.0 




0.0 


5 65 . 7 




0.0 


0.0 




0.0 


528 . 5 




0.0 


0.0 




0.0 


429.1 




0.0 


0.0 




0.0 


2 26.9 




0.0 


0.0 




0.0 
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TABLE XXIII 





DATA FILE 


DD- 


963 FOREMAST 
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20 


0.0 


0.0 


214.8 


0.0 


0.0 


0.0 


21 


0.0 


0.0 


468.9 


0.0 


0.0 


0.0 


22 


0.0 


0.0 


539.0 


0.0 


0.0 


0.0 


23 


0.0 


0.0 


539.0 


0.0 


0.0 


C . 0 


2A 


o.c 


0.0 


91.6 


0.0 


0.0 


0.0 


25 


0 . c 


0.0 


4 74 . 7 


0.0 


0.0 


0.0 


26 


0.0 


0 .0 


2234.2 


0.0 


c.o 


0.0 


27 


0 . 0 


0.0 


2296.2 


0.0 


0.0 


0.0 


28 


0.0 


0.0 


18 74.0 


0.0 


0 . 0 


0.0 


29 


0 .0 


0.0 


1874.0 


0.0 


0.0 


0.0 


30 


0.0 


0.0 


1459.2 


0.0 


0.0 


0.0 


21 


0 .0 


0.0 


1238.5 


0.0 


0 . 0 


0.0 


32 


o.c 


0.0 


1489.9 


0.0 


0.0 


0.0 


23 


0.0 


0.0 


211 . 7 


0.0 


0.0 


0.0 


34 


0 .0 


0 .0 


498.7 


0.0 


0.0 


c.o 


35 


o.c 


0.0 


606.8 


0.0 


0.0 


0.0 


36 


0 .0 


0.0 


498 . 7 


0.0 


0.0 


0.0 


27 


0 .0 


0.0 


725.9 


0.0 


0.0 


0.0 


38 


0.0 


0.0 


174.3 


0.0 


0.0 


0.0 


39 


0.0 


0.0 


725.9 


0.0 


0.0 


0 . 0 


AO 


o.c 


0.0 


174.3 


0.0 


0.0 


0.0 


A1 


o.c 


0.0 


345 . 5 


0.0 


0.0 


0.0 


A2 


0.0 


0.0 


1729.0 


0.0 


0.0 


0.0 


43 


0.0 


0.0 


1729.0 


0.0 


0.0 


0.0 


44 


0.0 


0.0 


14 2 5 . 6 


0.0 


0.0 


0.0 


45 


0 .0 


0.0 


1810.4 


0.0 


0.0 


0.0 


46 


0.0 


0.0 


1468 . 1 


0.0 


0.0 


0.0 


47 


0.0 


u.o 


1810.4 


0.0 


0.0 


G. 0 


46 


0.0 


0.0 


1468.1 


0.0 


0.0 


0.0 


49 


o.c 


0.0 


185.6 


0 . 0 


0.0 


0.0 


50 


0.0 


0.0 


185.6 


0.0 


c.o 


0.0 


51 


o.c 


0.0 


5 40.9 


0.0 


0.0 


0.0 


52 


0.0 


0.0 


965 . 3 


0.0 


0 . 0 


0.0 


53 


0.0 


0.0 


783.4 


0.0 


0.0 


0.0 


54 


0.0 


0.0 


464.2 


0.0 


0.0 


0.0 


55 


0.0 


0.0 


113.0 


0.0 


0.0 


0.0 


56 


0.0 


0.0 


1 12.0 


0.0 


0.0 


0.0 


57 


0.0 


0.0 


294 . 5 


0.0 


0.0 


0.0 


58 


0.0 


0.0 


436.2 


0.0 


c.o 


0.0 


59 


0.0 


0.0 


1425.6 


0.0 


0.0 


0.0 


61 


0.0 


0.0 


733 . 8 


0.0 


0.0 


0.0 


62 


0.0 


0.0 


304.6 


0.0 


0.0 


0.0 


63 


0.0 


0.0 


544 . C 


0.0 


0.0 


0.0 


64 


0.0 


0.0 


14C7.7 


0.0 


0.0 


0 . 0 


65 


0.0 


0.0 


125.0 


0.0 


0.0 


0.0 


66 


0.0 


0.0 


138 . 8 


0.0 


0.0 


0.0 


67 


0.0 


0.0 


100.3 


0.0 


0.0 


o.c 


68 


o.c 


0.0 


60 . 1 


0.0 


0.0 


0.0 


69 


0 .0 


0.0 


60 . 1 


0.0 


0.0 


0.0 


70 


0.0 


0.0 


463.2 


0.0 


0.0 


0.0 


71 


o.c 


0.0 


1713.1 


0.0 


0.0 


0.0 


72 


0.0 


0.0 


237.0 


0.0 


c.o 


C . 0 


73 


0.0 


0.0 


1081.2 


0.0 


0.0 


0.0 
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PLATFORM 3FLCW JPQ-S 



18 
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2C 

21 



32 .0 
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22 . C 
32*0 



SPC-S ANC FLATFCPM 



26 

27 

28 

29 

30 

31 
22 



590.0 
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890.0 
890. C 
890. C 
890.0 



SPG-60 ANC FUTFCRP 



NFMASS CAROS FORMAT I10 , F10.0 



90 



s 



A2 


540.0 


^3 


940. C 


44 


940 .0 


45 


940 .0 


46 


940. C 


47 


940.0 


48 


940.0 


59 


940 .0 


►55 A^O PLATFI 


22 


325.0 


23 


325 .0 


E4 


325. C 


iNNAS 


AhC OTH 


49 


25 .0 


5C 


25.0 


c t; 


12.0 


56 


12.0 


64 


220.0 


67 


30. C 


68 


26. C 


69 


26.0 


7C 


120 .0 


71 


972 . C 


73 


565. C 







56LCCK L DESIGN 


VARIABLE 
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AS R6Q*0 


FORMAT 


3 F 10.0 




S (AREA AND 


CIMENSICN 


VARIABLES 


GMI T 


IF N0VAR1=' 


0) 






S XA(I) . 


, XA<NOVARlj 


1 


1 


1 


1 




5 1 

13.75 


. 74 


7.625 


.37^ 


8.625 


.37^ 


7.625 


.375 


8.625 


.375 


12.50 


.50 


7.625 


.375 


7.625 


.27 5 


7.625 


.375 


23.00 


1.50 


9.50 


.50 


14.00 


.75 


10.00 


. 50 


7.625 


.375 


7. 625 


.375 


11 .50 


.50 


5.75 


.25 














i XALdj 


1 


1 


1 


1 


1 


1 




8. CO 


.50 


5.000 


.250 


5.000 


.250 


5 . 00 c 


. 250 


5.CC0 


.250 


8.00 


.25 


5.000 


.250 


5.0CC 


.250 


5.0CO 


.2 50 


15.00 


.50 


6.00 


.25 


10. OC 


.50 


6. CO 


.25 


5.000 


.250 


3.00 


.10 


3.00 


.25 


3. CO 


. 10 














$ XAUdj 






1 


1 


I 


1 




S 1 

22. CO 


1 

2. 50 


1 

15. CO 


1.50 


15.00 


1 

1.50 


15.00 


1.50 


15. CO 


1.50 


22.00 


2.50 


15.00 


1.50 


15. CO 


1.50 


15. CO 


1. 5C 


36.00 


2.50 


15.00 


1.50 


26.00 


2.00 


15.00 


1. 50 


15.00 


1.50 


12.00 


1.50 


22.00 


1.50 


12. CC 


1. 50 
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NCSPLC^C 


IN BLOCK 8> 












1 JN CIR LC 3L 




8U 








<S 1 

67 


1 


i 


1 

0 


6.6 








67 


3 


1 


-6. 0 


6.0 








46 


3 


1 


-3.0 


3.0 








48 


3 


1 


-3.0 


3.0 








22 


1 


1 


-6. 0 


6.0 








22 


3 


1 


-6.0 


6.0 








23 


1 


1 


-6. 0 


6.0 








23 
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APPENDIX 3 



PROGRAM ORGANIZATION 



A, DESCRIPTION 

The program organization is layed out in the following 
flow charts. The main driver program (SADXM) calls a 
subdriver (SADXSD), and the optimizer of the user's choice. 
All changes required for replacement of the optimizer are 
.made in SADXM, thus allowing for easy testing of several 
optimizers on the same problem. 

SADXSD may be called from the main for input, analysis, 
and output. Printed output may vary as the user requires. 

A complete listing of all subroutines and their functions is 
given in Table XXIV. 
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TABLE XXIV 

SUBROUTINE DIRECTORY 

DRIVER PR0S51H~?0R“aSING THE ABOVE 
SU3ROUTIHES. SADXH MAI BE COUPLED 
TO OPTIMIZER OF USER'S CHOICE. 


SADXSD 


SUBDRIVSR PROGRAM FOR COUPLING 
3ACX ROUTINES TO SADXM 


SADX01 


THIS ROUTINE READS AND PRINTS INPUT DATA 
AND ORGANIZES PS SDDO-DINAHIC STORAGE 
ALLOCATION 


SADX02 


BUILDS VECTORS JC AND IIK FOR FINITE ELEMENT 
STRUCTURAL ANALYSIS 


SADX03 


BUILDS THE 1 2x12 ELEMENT STIFFNESS MATRIX 


SADX05 


SUPERIMPOSES THE ELEMENT STIFFNESS MATRIX EK 
(OR ELEMENT MASS MATRIX EM) ON THE COMPACTED 
GLOBAL STIFFNESS MATRIX AK (OR THE GLOBAL 
COMPACTED MASS MATRIX AM) 


SADX06 


BUILDS GLOBAL LUMPED MASS MATRIX 


SADX07 


LU DECOMPOSES SY MMETP.IC, POSITIVE- DEFINITE, 
SPARCS MATRICES, THE UPPER TRIANGLE OF NHICH 
IS STORED IN MATRIX AK (OR AM) NITH LEADING 
ZEROES NOT STORED 


SADX08 


FORWARD AND BACK SUBSTITUTES TO YIELD A 
SOLUTION A SET OF LINEAR EQUATIONS 
(DECOMPOSED EY SADX07 DR EQUIVALENT) 


SADX09 


PRINTS ALL JOINT DISPLACEMENTS FOR EACH 
LOAD CONDITION OF A FINITE ELEM. STRUCTURE 


SADX1 1 


ROUTINE rO ORGANIZE ANALYSIS 


SADX15 


CALCULATES VALUES FOR ALL DFSIGN AND 
3E3AVIORIAL CONSTRAINTS AS DEFINED BY 
"SADX” PROGRAM 


SADX1 6 


CALCULATES STRESS IN TRUSS ELEMENT LNO 
UNDER LOAD CONDITION JJ 


SADX17 


PRINTS STRESSES AND/OR FORCES FOR 
TRUSS ELEMENTS 


SADX19 


ADDS ELEMENT MASS MATRIX A A OF ELFMSNT 
LNO TO GLOBAL MASS MATRIX AM TO BUILD 
THE LUMPED MASS MATRIX 


SADX23 


CALCULATES WEIGHT OP IRUSS/PRAME STRUCTURE 
OR CALCULATE WEIGHT DF INDIVIDUAL MEMBERS 


SADX36 


CALCULATES ( XS IG-T’''AM* XEI G) FOR 
GRADIENT CALCULATIONS IN FREQUENCY 
CONSTRAINTS 


SADX37 


CALCULATES EIGENVALUE GRADIENT INFORMATION 
IN FINITE ELEMENT STRUCTURAL ANALYSIS AND 
DESIGN 
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SADX46 

SADX47 

SADX49 

SADX50 

SADX53 

SADX62 

SADX71 

SADX72 

SADX78 

SADX80 

SADX81 

SADX82 

SADX83 

SADX84 

SADX85 

SADX86 

SADX87 



HSADS INPUT INFORMATION FOR TRUSS 
ELEMENTS 



TRANSFORMS THE ELEMENT STIFFNESS MATRIX 
EK (OR ELEMENT HASS MATRIX EM) FROM LOCAL 
TO GLOBAL COORDINATES 



SOLVES REAL EIGENVALUE PROBLEMS USING 
THE SUES? ACE ITERATION METHOD 



BUILDS INITIAL SET OF BASIS VECTORS 
FOR SIGENSOLUTION BX REDUCED BASIS 
METHOD 

PRINTS MEMBER INFORMATION FOR TRUSS 
ELEMENTS 

PRINTS NEIG EIGENVALUES STORED IN 
EIGVAL, AND THEIR CORRESPONDING 
EIGENVECTORS STORED IN XEIG 

PRINTS G VECTOR OF CONSTRAINTS 

READS IN FRAME ELEMENT INPUT DATA 

BUILDS 3x3 TRANSFORMATION ARRAY TRFORM 
FOR TRANSFORMING FROM LOCAL TO GLOBAL 
COORDINATES 

CALLS SADX03 TO BUILD THE ELEMENT STIFFNESS 
MATRIX; CALLS SADX47 TO TRANSFORM THE MATRIX 
AND CALLS SADX05 TO BUILD THE REDUCED GLOBAL 
STIFFNESS MATRIX 

CALCULATES STRESS IN FRAME ELEMENT LNO UNDER 
LOAD CONDITION JJ (IF JJ=0 STRESSES CALCUL- 
ATED FOR ALL LOAD CONDITIONS) 

READS INPUT DATA FOR FRAME ELEMENTS WITH 
SPECIFIED SECTION TYPES 

PRINTS STRESSES AND/OR FORCES FOR FRAME 
ELEMENTS 

PRINTS MEMBER INFORMATION FOR FRAME ELEMENTS 

CALCULATES SECTION PROPERTIES FOR FRAME 
ELEMENTS OF A SECTION TYPE GIVEN BY LSECT 

CALLS EITHER SADX06 TO BUILD THE LUMPED MASS 
MATRIX OR BUILDS THE CONSISTENT HASS MATRIX 
BY CALLING SADX87 TO BUILD THE ELEMENT MASS 
MATRIX, SA0X78 TO BUILD THE TRANSFORMATION 
MATRIX, SADX47 TO TRANSFORM THE ELEMENT MASS 
MATRIX, AND SADX88 TO ASSEMBLE THE COMPACTED 
GLOBAL MASS MATRIX 

CALLS BUILDS THE ELEMENT CONSISTENT MASS 
MATRIX 
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•“j; 



SADX88 

SADX89 

SADX90 

SADX91 

SETISE 

GSTIME 



CONVESTS UNFORMATTED DATA TO FORMATTED DATA 
IN FIELDS OP 10, EACH RIGHT JUSTIFIED AND 
ACCEPTS SCOaHENT CARDS IN DATA 

SOLVES EIGENVALUE PROBLEM |A-ALAMBD*Bi 7=0 

SOLVES EIGENVALUE PROBLEM 

SOLVES EIGENVALUE PROBLEM 

STARTS EXECUTION TIMER 
(NONiaSL LIBRARY) 

STOPS EXECUTION TIMER 
(NONIMSL LIBRARY) 
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TABLE XXV 



PROGRAM BLOCK DIAGRAM 




SADX85 



ADX03 

ADX78 

ADX47 

ADX05 

ADX06- 

ADX87 

ADX78 

ADX47 

ADX05 



ADX08 

ADX89> 



ADX78 

ADX03 

ADX78 

ADX03 



ADX73 
ADX03 
A 0X73 
ADX03 



SADX 19 



SADX90 

3ADX91 
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